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c§ ■ ABSTRACT 

> : 

^ \ We determine the primordial helium mass fraction Yp using 93 spectra of 86 

\ low-metallicity extragalactic H ll regions. This sample constitutes the largest 

and most homogeneous high-quality data sets in existence for the determination 
of Yp. For comparison and to improve the statistics in our investigation of sys- 
tematic effects affecting the Yp determination, we have also considered a sample 
of 271 low-metallicity H ll regions selected from the Data Release 5 of the Sloan 
Digital Sky Survey. Although this larger sample shows more scatter, it gives 
results that are consistent at the 2a level with our original sample. We have 
considered known systematic effects which may affect the ^He abundance deter- 
mination. They include different sets of He I line emissivities and reddening laws, 
collisional and fluorescent enhancements of He I recombination lines, underlying 
He I stellar absorption lines, collisional excitation of hydrogen lines, temperature 
and ionization structure of the H ll region, and deviation of He I and H emission 
line intensities from case B. However, the most likely value of Yp depends on the 
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adopted set of He I line emissivities. Using Monte Carlo methods to solve simul- 
taneously the above systematic effects we find a primordial heli um mass fraction 



Yp = 0.2472 ± 0.0012 when using the He I emissivities from iBeniamin et al 



(119991, 12OO2I) and 0.2516 ± 0.0011 when using those from IPorter et al. I (120051 ). 



The first value agrees well with the value given by Standard Big Bang Nucle- 
osynthesis (SI 3BN) theory, while the value obtained with likely more accurate 
emissivities of iPorter et al. I (120051 ) is higher at the 2a level. This latter value, if 
confirmed, would imply slight deviations from SBBN. 



Subject headings: galaxies: abundances 
— H II regions — ISM: abundances 



galaxies: irregular — galaxies: ISM 



INTRODUCTION 



In the standard theory of big bang nucleosynthesis (SBBN), given the number of light 
neutrino species, the abundances of these light elements depend only on one cosmological 
parameter, the baryon-to-photon number ratio 77, which in turn is directly related to the 
baryon density parameter the present ratio of the baryon mass density to the critical 
density of the Universe. This means that accurate measurements of the primordial abun- 
dances of each of the four light elements can provide, in principle, a direct measurement of 
the baryonic density. As 77 is a very small quantity, it is convenient to define the parameter 
TjiQ = 10^° rj. Then rjio is related to Qb by the expression rjiQ = 274 Qh h'^ , where h = ff n/100 
km s~^ Mpc~^ and Hq is the present value of the Hubble parameter (jSteigman II2005I ). 



Because of the strong dependence of its abundance on i], deuterium has become the bary- 
ometer of choice ever since accurate measurement of D/H in high-redshift low-metallicity 
QSO Lya absorption systems have become possible. Although the data is still scarce - 
there are only six absorption syst e ms for which such a D / H measurement has been car 



ried out (Buries fc Tvtler 



Kirkman et al. 



1998a 



b: O'Meara et al. 2001, 2006: Pettini fc Bowen 2001 



2OO3I : Crighton et al. 2004) - and the scatter remains large, the measure 



ments appear to converge to a mean primordial value D/H ~ 2.4 ± 0.4 x 10~^, which corre- 
sponds to Qbh"^ ~ 0.023 ± 0.002. This estimate of flbh"^ is in excellent agreement with the val- 
ues of 0.022 - 0.023 obtained from recent studies of the fluctuations of the cosmic microwave 



and WMAP (ISpergel et al. 112003 



background ( CMB) bv DASI (IPrvke et al. II2002I) . BOOMERANG (INetterfield et al. II2OO2I ) 



200' 



derived from large scale structure (ITegmark et al. II2004I ). 



^). It is also in goo d agreement with the value of ilbh"^ 



While a single good baryometer like D is sufficient to derive the baryonic mass den- 
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sity from BBN, accurate measurements of the primordial abundances of at least two dif- 
ferent relic elements are required to check the consistency of SBBN. Among the remaining 
relic elements, ^He can only be observed in the solar system and in the Galaxy, both of 
which have undergo ne significant chemical evolution, making it difficult to derive its pri- 
mordial abundance (IBania et al. 1 12002| ) . The derivation of the primordial ^Li abundance 
in metal-poor halo stars in the Galaxy is also beset by difficulties such as the uncertain 
stellar tem perature scale and the temp erature structures of the atmospheres of these very 
cool stars (jCharbonnel fc Primas II2005I ). The primordial abundance of the remaining relic 
element, ^He (hereafter He), can in principle be derived accurately from observations of the 
helium and hydrogen emission lines from low-metallicity blue compact dwarf (BCD) galax- 
ies which have undergone little chemical evolution. Several groups have used this technique 
to derive the primordial He niass f r action Y^, with somewhat diff erent results. Our group 



fllzotov. Thuan fc Lipovetskv 1119971 : llzotov fc Thuan Ill998bl . 120041 ) has obtained Yp = 0.244 



± 0.002. This va lue is consistent at the two-sigma lev el with the prediction of SBBN, Yp 



0.2482 ± 0.0007 flWalker et al. 



1991 



Steigman 



± 0.0008 found by WMAP fISpergel et al. 1120061 ). On the other hand, the Fp predi c ted b y 



20051 ) ■ adopting the value of Qbh"^ = 0.0223 



SBBN is not consistent with the lower values o f 0.234 ± 0.002 derived bv lOlive et al. I ( 



1997), 



of 0.2 384 ± 0.0025 by iPeimbert et al. I (120021 ) and of 0.2391 ± 0.0020 by iLuridiana et al 
(120031 ). These lower estimates of Yp have led some authors to argue that the discrepancy 



between the flhS derived from Yp and from D/H is a sign that SBBN is not valid (jChen et al 



200ll : iKneller et al. II2OOII : iBarger et al. ll2003al Jbl: ICvburt et al. II2005I ) and that new physics 
are required. Although He is not a sensitive baryometer [Yp depends only logarithmically on 
the baryon density), its primordial abundance depends much more sensitively on the expan- 
sion rate of the Universe and on a possible asymmetry between the numbers of neutrinos and 
anti-neutrinos in the early universe than D, ^He or ''Li. This, because of two reasons: 1) a 
faster expansion would leave less time for neutrons to convert into protons, and the resulting 
higher neutron abundance would result in a higher Yp, and 2) Yp depends sensitively on the 
neutron to proton ratio, which depends in turn on the numbers of electron neutrinos and 
anti-neutrinos. In that sense. He is both a chronometer and/or leptometer which is very 
sensitive to any small deviation from SBB N, and hence to new physics, much more so than 
the other three primordial light elements (jSteigman II2006I ). 



However, to detect small deviations from SBBN and make cosmological inferences, Yp 
has to be determined to a level of accuracy of less than one percent. This is not an easy task. 
While it is relatively straightforward to derive the helium abundance in an H ll region with 
an accuracy of 10 percent if the spectrum is adequate, gaining one order of magnitude in 
the precision requires many conditions to be met. First, the observational data has to be of 



excellent quality. This has been the concern of our group (llzotov. Thuan &: Lipovetsky Ill994 
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19971 : llzotov fc Thuan Ill998bl . l2004j ). We have spent the last decade obtaining high signal-to- 
noise spectroscopic data of low-metallicity extragal actic H ll regions, and o ur sample includes 
now a total of 86 H ll regions in 77 galaxies (see llzotov fc Thuan 1 120041 ) . This constitutes 
by far the largest sample of high-quality data reduced in a homogeneous way to investigate 
the problem of the primordial h elium abundance. To put things in perspe ctive, the sample 



used in the pioneering work of iPeimbert fc Torres-Peimbert I (Il974j . (19761) co mprised only 
5 obje cts, with considerably larger ob servational errors. Later, the sample of iPagel et al. 
( I1992I ) included 36 objects and that of lOlive fc Steigman I (119951 ) 49 objects. 



With such a large observational sample at hand, it is now the general consensus that 
the accuracy of the determination of the primordial He abundance is limited presently, not 
so much by statistical uncertainties, but by our ability to account for systematic errors and 
biases. There are many known effects we need to correct for to transform the observed 
He I line intensities into a He abundance. Neglecting or misestimating them may lead to 
systematic errors in the Yp determination that are larger than the statistical errors. These 
effects are: (1) reddening, (2) underlying stellar absorption in the He I lines, (3) collisional 
excitation of the He I lines which make their intensities deviate from their recombination 
values, (4) fluorescence of the He I lines which also make their intensities deviate from their 
recombination values, (5) collisional excitation of the hydrogen lines (hydrogen enters because 
the helium abundance is calculated relative to that of hydrogen), (6) possible departures 
from case B in the emissivities of H and He I lines, (7) the temperature structure of the 
H II region and (8) its ionization structure. All these corrections are at a level of a few 
percent except for effect (3) that can be much higher, exceeding 10% in the case of the He 
I A5876 emission line in hot and dense H ll regions. All effects, apart from (8), influence 
each other in a complicated way. At the present time, we are far from a situation where 
all these corrections can be determined unambiguously from observational constraints and 
theory. Given the complexity of H ll regions - they have very perturbed morphologies, 
possess complex velocity fields and a dust distribution difficult to quantify - the situation is 
not likely to be improved considerably in the near future. Because of this, even if one were 
able to solve in an exact manner the radiation transfer in the case of a simple geometry, 
there is still a long way between a simple model and a real H ll region. 

Effects (1), (3), (4) are relatively easy to account for empirically from the available 
spectroscopic information an d have been already considered in the pioneering study of 



Peimbert fc Torres-Peimber 



recently (IRayo et al. 



1982 



(Il974r) . Corrections (2) and (5) have been disc ussed more 



Kunth fc Sargent Ill983l : IPavidson fc Kinman Ill985[ ). Effect (7) 



has been the subject of many papers by Peimbert and his coworkers. Elemental abundance 
determinat ions are affected by bi ases in the presence of te mperature inhomogeneities in the H 



II regions. iPeimbert I (119671 ) and iPeimbert et al. I (120021 ) have devised an empirical method 



- 5 - 



to correct for this bias. However, this method rehes on several assumptions that are not 
necessarily valid (e.g. supposing that th e temperature fluctua tion parameter is the same 



in the high and low ionization regions, iPeimbert et al. II2002I ). Furthermore, the source 



of temperatu re inhomogeneities at the level invoked by Peimbert and colleagues remains 



unidentifled (IStasihska fc Schaerer I Il999l : iBinette fc Luridiana 



they are supposed to result from abundance inhomogeneities (ITsamis fc Pequignot 1 12005 



2000 



Mathis 1995), unless 



Stasihska et al., in preparation). This is also why the computation of temperature cor 



rectio ns from photoionization models to derive He/H (as proposed by ISauer fc Jedamzik 



20021 ), is not necessa r ily relevant to real objects. To check for possible temperature fluctua- 



tions, 



Guseva et al. 



(I2OO6I . I2OO7I ) have compared the temperature in the 0^+ zone derived 
with the [O ill] A4363 line, with that in the H"*" zone derived with the Balmer and Paschen 
jumps in a large sample of hot low-metallicity H ll regions used in the determination of Yp. 
They found that the two temperatures do not differ statistically, so that any temperature 
difference must be less than ~ 5%. Effect 



ical method (Peimbert & Torres- Peimbert 


1974). and has been 


ater the subiect of many 


papers based on photoionization modeling ( 


Stasim 


[a 


1980: Pefial 


1986:IVieeas et al. 


2000: 


Ballantvne et al. 2000: Sauer & Jedamzik 


2002: 


Gruenwald et al. 2002). Effect (6) has 


been discussed for H ll regions bv Cota & Ferland 


- ( 


1988) and Hummer & Storev f 


1992). 



but has never been considered in the various estimates of Yp thus far, case B having been 
assumed by all authors. 

In the present determination of the primordial He abundance, we take into account all 
eight effects discussed above. We flrst use an empirical method to self-consistently account 
for effects 1, 2, 3, 4, 5 and 7 in the deri vation of the abundance of He"*". This method has 
been detailed in llzotov fc Thuan I (120041 ). where it has been applied to a small sample of 7 
objects. It was shown that taking all these effects into account in a system atic way can lead 



to heh um abundances that are signiflcantly higher than those obtained by iPeimbert et al. 



(I2OO2I ) for the same objects (by up to 3%). Here, we apply the sa me method, not only to 
our basic original sample of 93 H ll regions (llzotov &: Thuan II2004I ) that we have assembled 
and observed ourselves (hereafter the HeBCD sample), but also to a different and larger 
sample of 271 extragalactic H ll regions selected from the Data Release 5 (DR5) of the Sloan 
Digital Sky Survey (hereafter the SDSS sample). We then considered effects 6 and 8 on the 
derivation of the abundance of He in each H ll region and of the primordial He abundance. 

We describe the two HeBCD and SDSS samples in ^ In ^ we discuss the method 
used to derive He abundances in individual objects and the primordial He abundance from 
each sample. In §11 we discuss the systematic effects considered. Our new best values for Yp 
and the linear regression slope dY/dZ are presented in ||5l The cosmological implications of 
our new results are discussed in J|6l We summarize our conclusions in ^ 
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2. THE SAMPLE 

2.1. The HeBCD and SDSS subsamples 

We determine Yp and dY/dZ independently for two different samples. The HeBCD 
sample, is composed of 93 different observations of 86 H ll regions in 77 galaxies. The ma- 
jority of these galaxies are low-metallicity BCD galaxies. This sample is essentially the same 



as the one described in llzotov fc Thuan I (j2004j ). except that two close pairs of H ll regions 
in the spiral galaxy MlOlIl have been added. These four H ll regions were observed with 
the 6.5m MMT on the night of 2004 February 20. Observations were made with the Blue 
Channel of the MMT spectrograph. We used a 2"x300" slit and a 300 grooves/mm grating 
in first order. To remove the second-order contamination we use the L-38 blocking filter. The 
above instrumental set-up gave a spatial scale along the slit of 0'.'6 pixel"^, a scale perpen- 
dicular to the slit of 1.96A pixel"^, a spectral range 3600 - 7500A and a spectral resolution 
of 6A (FWHM). The Kitt Peak IRS spectroscopic standard star HZ 44 was observed for flux 
calibration. Spectra of He-Ar comparison arcs were obtained after each o bservation to cali- 



brate the wavelength scale. The details of data reduction are described in iThuan fc Izotov 



( 120051 ) and they are the same as those for other H ll regions from the HeBCD sample. The 



two-dimensional spectra were bias subtracted and fiat-field corrected using We then 

use the IRAF software routines IDENTIFY, REIDENTIFY, FITCOORD, TRANSFORM to 
perform wavelength calibration and correct for distortion and tilt for each frame. Night sky 
subtraction was performed using the routine BACKGROUND. The level of night sky emis- 
sion was determined from the closest regions to the galaxy that are free of galaxian stellar 
and nebular line emission, as well as of emission from foreground and background sources. 
One-dimensional spectra were then extracted from each two-dimensional frame using the 
APALL routine. Before extraction, distinct two-dimensional spectra of the same H ll region 
were carefully aligned using the spatial locations of the brightest part in each spectrum, so 
that spectra were extracted at the same positions in all subexposures. For all objects, we 
extracted the brightest part of the H ll region. We have summed the individual spectra from 
each subexposure after manual removal of the cosmic rays hits. The spectra obtained from 
each subexposure were also checked for cosmic rays hits at the location of strong emission 



^ The coordinates of the two H ii region pairs are: for the pair MlOl No.l + No. 2, R.A.(J2000.0) 
14''04™29!5, Dec (J2000.0) = +54°23'47" and for the pair MlOl No.3 + No.4, R.A.(J2000.0) = 14h03"^01!2, 
Dec (J2000.0) = +54°14'29". 

^IRAF is distributed by National Optical Astronomical Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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lines, but none was found. 

Particular attention was paid to the derivation of the sensitivity curve. It was obtained 
by fitting with a high-order polynomial the observed spectral energy distribution of the 
standard star HZ 44. Because the spectrum of HZ 44 has only a small number of a relatively 
wea k absorptio n features, its spectral energy distribution is known with good accuracy of < 



1% (lOke lllQQOl ). Moreover, the response function of the CCD detector is smooth, so we could 



derive a sensitivity curve with an accuracy better than 1% over the whole optical range. 

We show in Tables [1] and [2] (available only in electronic form), the emission line fluxes and 
the equivalent widths for these H ll regions, along with those of the remaining H ll regions 
in the HeBCD sample. They were measured using the IRAF SPLOT routine. The line flux 
errors listed include statistical errors derived with SPLOT from non-fiux calibrated spectra, 
in addition to errors of 1% of the line fluxes introduced in the standard star absolute flux 
calibration. The line flux errors will be later propagated into the calculation of abundance 
errors. 

The number of data points (93) is larger than the number of H ll regions (86) be- 
cause several H ll regions have independent observations from different telescopes. We 
treat these independent observations as separate data points in our least-square fitting. 
In assembling the above sample, we have taken care not to include objects that, for one 
reason or another, are not appropr i ate fo r He abundance deter n ainatio n, as described in 
Izotov. Thuan fc Lipovetsky I (11994 119971 ) and llzotov &: Thuan I (ll998al Jbl). In particular. 



the NW component of I Zw 18 has not been included because its He I A5876 emission line is 
strongly affected by the Galactic interstellar Na I A5889,5895 absorption line. 

The SDSS sample, which we will use as a comparison sample, is composed of 271 



low-metallicity H ll regions selected from the SDSS DR5. The SDSS (lYork et al. 1 120001 ) 
offers a gigantic data base of galaxies with well-defined selection criteria and observed in 
a homogeneous way. In addition, the spectral resolution is much better than that of most 
previous data bases on emission-line galaxies including all the spectra in the HeBCD sample. 
It is possible to extract from the SDSS data base a sample of emission-line galaxies with 
well-defined criteria. Despite the lower quality of its data, this large comparison sample, 
observed and reduced in a different way, will allow us to check for possible systematic shifts 
in the linear regressions introduced by different data sets. Its large size will also allow to 
check for other systematic effects. 

The SDSS DR5 provides spectra of some 800 000 galaxies, quasars and stars selected 
over a sky area of 4783 square degrees, in the wavelength range ~ 3800 - 9200A, along with 
tables of measured parameters from these data. From this data base accessible from the 
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SDSS web page (http://www.sdss.org/dr5), we have extracted flux-cahbrated spectra of a 
total of 271 H II regions which satisfy the following selection criteria: 1) the [O ill] A4363 
is detected at the level of > 2o" above the noise, allowing a direct heavy element abundance 
determination by the Te method, as was done for the H ll regions in the HeBCD sample; 2) 
the equivalent width of the H/3 emission line EW(H/5) is > 50A. This ensures that the effect 
of underlying He I stellar absorption on the derived He abundances is small; 3) the observed 
flux of the H/3 emission line F{B.(3) is > 10"^*^ erg s~^ cm~^. This condition ensures that the 
SDSS sample contains only bright H ll regions with strong emission lines, the fluxes of which 
can be measured accurately. Since the emission line [O ll] A3727 in the SDSS spectra of 
galaxies with z < 0.02 is out of the observed spectral range, we have estimated its flux from 
the [O ll] A7320, 7330 emission line fluxes. The measurements of the emission line fluxes 
and the determination of the element abundances were done following the same procedures 
as for our basic HeBCD sample. 

The observed emission line fluxes F{X)/F{}if3) and their equivalent widths EW(A) for 
all the objects in the two samples are shown respectively in Tables [1] and [2l Because of their 
large sizes, these tables are available only in electronic form. 



2.2. Statistical errors 

Emission line fluxes of all spectra, in both the HeBCD and SDSS samples, were measured 
using Gaussian fitting with the IRAF routine SPLOT. The statistical errors of emission lines 
in the spectra of the HeBCD sample are calculated using the photon statistics of the lines in 
the non-flux calibrated spectra. As for the SDSS sample, each spectrum in the SDSS data 
base is supplemented by a file with the flux error in each pixel which can be used to calculate 
the errors of the emission lines in each spectra. The mean statistical errors are ~ 4% for the 
brightest He I A5876 emission line, ~ 10% for the He I A4471 and A6678 emission lines, and 
~ 13% for the He I A 7065 emission line, for the combined HeBCD + SDSS sample. As the 
mean statistical error of the weighted mean He abundance is primarily determined by the 
statistical error of the intensity of the brightest He I A5876 emission line, we expect it to be 
~ 3-4%, which is the case (see Table [3]). 

To check the reliability of error determinations for the emission lines in the spectra of 
our sample, we have considered the flux ratios of the two [O l] A6363 and [O l] A6300 lines. 
These flux ratios are plotted against oxygen abundance in Fig. [T]where objects in the HeBCD 
sample are shown by filled circles, and those in the SDSS sample by open circles. The solid 
line shows the mean value of the flux ratio for the combined HeBCD + SDSS sample, and the 
dashed lines represent la dispersions, ~ 17% of the mean value. Since the theoretical value 
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of the [O l] A6363/A6300 ratio is constant (its value is ~ 1/3) and independent of physical 
conditions in the H II region, the dispersion of the points around the mean in Fig. [1] gives a 
good representation of the observational uncertainties. Typically, the flux of [O l] A6300 is 
comparable to the flux of He I A4471 and A6678 emission lines, with statistical errors of ~ 
10%. The [O l] A6363 is ~ 3 times weaker than the [O l] A6300 line, so its error is ~ 17%. 
This is just the la scatter found in Fig. [T]for the [O l] line flux ratio. We conclude that our 
estimates of the statistical errors for weak emission lines are reliable. 

The statistical errors in the emission lines intensities are then propagated in the calcu- 
lation of errors in electron temperatures, electron number densities and element abundances. 



3. THE METHOD 



3.1. Linear regressions 



As in our previous work (see llzotov fc Thuan 1 12004| . and references therein), we de- 
termine the primordial He mass fraction Yp by fitting t he data points in the Y - 0/B. and 
y - N/H planes with lir iear regression lines of the form (jPeimbert fc Torres-Peimbert Ill974j . 



19761 : IPagel et al. Ill992[ ) 



where 



Y = Y^ + 



Y = Yp + 



Y 



dY 



d(0/H) 

dY 
d(N/H) 

Ay{l-Z) 



(0/H), 



(N/H), 



(1) 



(2) 



(3) 



1 + Ay 

is the He mass fraction, Z is the heavy element mass fraction, y = {y^ + y^^) x/C-F(He"'"+He^"'") 
is the He abundance, = He+/H+ and y"^^ = He^+/H+ are respectively the abundances 
of singly and doubly ionized He, and JCF(He^-|-He^^) is the ionization correction factor 
for He. We have assumed Z = 18.2(0/H) in Eq. [3] which holds for a metallicity Z=0.001 



(IMaeder Ill992h . 



The assumption of a linear dependence of Y on 0/H and N/H appears to be reasonable 
as there are no evident non-linear tre nds in the distributions of the data points in the Y vs 
0/H and Y vs N/H diagrams (e.g., llzotov &: Thuan 1 120041 ) . The linear regressions (Eqs. 



[T] and [2]) imply that the initial mass function (IMF) averaged stellar yields for different 
elements do not depend on metallicity. This is the case i f the stellar IMF i s independent 



of metallicity. It has been suggested in the past (e.g., iBond et al. I Il983l ) that, at low 
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metallicities, the IMF may be top-heavy, i.e. there are relatively more massive stars as 
compared to lower mass stars than at high metallicities. If this is the case, then the IMF- 
averaged yields would be significantly different for low-metallicity stars as compared to those 
of more metal-enriched stars, result ing in a non-linear relationship between Y and 0/H 
or N/H (jSalvaterra &: Ferrara 1 120031 ) . However, until now, there has not been persuasive 
evidence for a metallicity dependence of the IMF. Furthermore, the properties of extremely 
metal-deficient stars remain poorly known, excluding quantitative estimates of possible non- 
linear effects in the Y - O/H and Y - N/H relations. Therefore, in the following analysis, 
we will continue to use linear regressions (Eqs. [1] and [2]) to fit the data. 



The slopes of the Y - O/H and Y - N/H linear regressions can be written as: 



dY 



d(0/H) 



dY 



dY 



dY 



10.5- 



dY 



564- 



.dY 



(4) 



(5) 



d(N/H) "■"'dN ""'dZ 

where O, N and Z are respectively th e mass fractioii s of oxygen, nitrogen and heavy elements. 
We have assumed that O = 0.66Z ( iMaeder Ill992h which holds for a metallicity Z=0.001, 
an IMF s lope a;=1.35 (where x is defined by dA^/d(logM) oc M^) and log(N/H) - log(0/H) 
= -1.55 jThuan et al. Ill995h . 



To derive the paramet ers of the linear regressions, we use the maximum-likelihood 
method (jPress et al. I Il992l ) which takes into account the errors in Y, O/H and N/H for 
each object. 



3.2. He+ emissivities 



The derived He"*" abundance depends on t he adopted He I l i ne em i ssivit ies. We con 



sider tw o sets of He I emissivitie s: the old ones by 



used by llzotov fc Thuan I (120041 ) Benjamin et al. 



and fl uorescent enhancements] and the new ones by iPorter et al. I ( 120051 ) and iBauman et al. 



Beniamin et al. I (119991 . |2002| ) which were 



(2002) take into account both coUisional 



which have been computed usi i ig ini prove d radiative and coUision al data. Follow- 



ing 



Izotov. Thuan &: Lipovetsky I (Il994j . 119971 ) and llzotov &: Thuan I ( Il998bl ) , we use the five 



strongest He I A3889, A447 1 , A587 6, A6678 and A7065 emission lines to derive Ne{B.e'^) and 



r(A3889). iBauman et al. I (120051 ) have estimated the accuracy of new emissivities in the 
low-density limit and found that accuracy is better than 1% for He I A4471, A5876, A6678 
and A7065 emission lines, but it is not as good for the He I A3889 emission line. The He I 
A3889 and A7065 lines play an important role because they are particularly sensitive to both 



- 11 - 



quantities. Since the He I A3889 line is blended with the H8 A 3889 line, w e have subtracted 
the latter, assuming its intensity to be equal to 0.107 /(H/5) (lAUer 1 1 19841 ). In our spectra, 
other He I emission lines are seen, most often He I A3820, A4387, A4026, A4921, A7281. How- 
ever, we do not attempt to use these lines for He abundance determination because they are 
much weaker as compared to the five brightest lines, and hence have larger uncertainties. 



We have used the simple and conv e nient fits provided by iBeniamin et al. I (120021) to 
calculate y~^s with the iBenjamin et al. I (119991 ) emissivities. For the iPorter et al. I (120051 ) 
emissivities, we have assumed the functional dependence of these emissivities on and 
r(He I A3889), which r esults from the co l lision al and fluorescence enhancements, to be the 
same as the one for the iBenjamin et al. I (12002! ) emissivities. Eqs. |6] - [10] give the linear fits 
we have adopted for the ratios of the two sets of emissivities in the temperature range Te= 
10^ - 2x10^* and an electron number density = 100 cm~^, for each of the five lines: 



?/+(A3889) = y+(A3889)/(1.079 - 0.052 x t), 

t/+(A4471) = ?/+(A4471)/(1.020 - 0.026 x t), 
7/J(A5876) = yj^(A5876)/(0.956 + 0.011 x t), 
?/+(A6678) = y+(A6678)/(0.938 + 0.028 x t), 
?/J(A7065) = i/J(A7065)/(1.051 - 0.040 x t), 



(6) 

(7) 
(8) 
(9) 
(10) 



where Vp and Up are respecti vely He^ abundances calculated with iBenjamin et al. I (119991 ) 



(B) and iPorter et al. I (120051 ) (P) emissivities and t = 10 ^ Te. We found th at, with 



Porter et al. i (l2005i ) emissivities. He"*" abundances are higher than those derived with iBenjamin et al. 
(1l999h emissivities by about 0-2 percent for the A4471 line, and by about 5-6 percent for 
the A5876 line, in the range t = 1 - 2. 



3.3. A Monte Carlo algorithm for determining the best value of 



In addition to the emissivities, the derived abundances depend also on a num- 
ber of other parameters: the fraction A/(Hq;)//(Hq;) of the Ha emission line flux due 
to collisional excitation, the electron number density Ne{B.e'^), the electron temperature 
Te(He+), the equivalent widths EWa6s(A3889), EWafes(A4471), EWab.(A5876), EWafe,(A6678) 
and EWafes(A7065) of He I stellar absorption lines, and the optical depth r(A3889) of the He 
I A3889 emission hne. To determine the best value of ?/,t.^ (defined in Eq. [T^D, we use the 



Monte Carlo procedure described in ilzotov &: Thuan I (120041 ). randomly varying each of the 
above parameters within a specified range. First, we take into account collisional excitation 
effects for hydrogen. The value of the fraction of the Ha flux due to collisional excitation 
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is randomly generated 100 times within an adopted range. The fraction of the H/3 emission 
hne flux due to the coUisional excitation is adopted to be three times less than that of the 
Ha flux. For each generated fraction, the fluxes of the Ha and H/? lines due to the collisional 
excitation are subtracted from the total observed fluxes and then all emission line fluxes are 
corrected for interstellar extinction and element abundances are calculated. 

To calculate we vary simultaneously and randomly N^iiie'^), Te(He"'") and r(A3889) 
within their respective adopted ranges. We make a total of 10^ such realizations for every 
H II region, for a given fraction of the Ha emission line flux due to collisional excitation. 
Thus, the total number of Monte Carlo realizations we have performed for each H ll region 
is 100 X 10^ = 10'^. As for the He I underlying stellar absorption, we assume fixed values for 
EW„6,(A4471), chosen to be between and O.SA, and for the EWab.(A3889)/EWa6,(A4471), 
EW„6,(A5876)/EW„,,(A4471), EW,b,(A6678)/EW,,,(A4471), and EW,b,(A7065)/EW,,,(A4471) 
ratios. 

For each H ll region, we find the best solution for ?/+^ in the multi-parameter space 
defined above by minimizing the quantity 

where y^' is the He"*" abundance derived from the flux of the He I emission line labeled i, 
and cr{yf ) is the statistical error of y^. The quantity ?/+^ is the weighted mean of the He"*" 
abundance as given by the equation 

We use all five He I emission lines to calculate n = 5), but only three lines. He I 

A4471, A5876 and A6678 to compute {k = 3). This is because the fluxes of He I A3889 
and A7065 emission lines are more uncertain as compared to the other three He I emission 
lines. 

An estimate of the IcJsys systematic error of ?/^^ can be obtained from Ax^. This 
quantity depends on the number of degrees of freedom in the problem which is the difference 
between the number of observational constraints, equal to 5 (the five He I emission line 
fluxes), and the number of the free parameters, equal to 4 [the fraction A/(Ha)/J(Ha) due 
to collisional excitation of the Ha emission line, Te(He+), A^pf He"*") and r(A3889) ]. Since the 



number of degrees of freedom is 5 - 4 = 1, then Ax^ = 1 fjPress et al. Ill992l ). Thus the 



systematic error osys is the la dispersion of the computed values of ?/+^ for solutions with 
between Xmm and Xmm + 1- The total error for ?/+^ is then given by a^^^ = a^^^ + a^^,. 
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Additionally, in those cases when the nebular He ll A4686 emission line was detected, 
we have added the abundance of doubly ionized helium y^^ = He^"'"/H'^ to Although 
the He^^ zone is hotter than the He"*" zone, we have adopted Te(He^"^) = Te(He"'"). The last 
assumption has only a minor effect on the y value, because is small (< 3% of y^) in all 
cases. 



3.4. A basic set of parameters for Monte Carlo calculations 

In order to study the dependence of Yp on the various parameters characterizing different 
systematic effects, we define a reference set of parameters which we will call the "basic set" . 
We wish to see how, by changing different parameters, the y values change, as compared to 
those corresponding to the basic set. The basic set of parameters is defined in the following 
way : 1) the electron temperature of the He"*" zone is varied in the range Te(He"'") = (0.95 

- 1.0)xTe(O III); 2) oxygen and nitrogen abundances are calculated adopting an electron 
temperature equal to T^iO ill); 3) Nf,{)^e^) and r(A3889) vary respectively in the ranges 10 

- 450 cm~'^ and 0-5; 4) the fraction of Ha emission due to collisional excitation is varied 
in the range 0% - 5%; 5) the equivalent width of the He I A4471 absorption line is fixed 
to EWa;,s(A4471) = 0.4A; 6) the equivalent widths of the other absorption lines are fixed 
according to the ratios EWab,(A3889)/EWai„(A4471) = 1.0, EWafe,(A5876)/EW„b,(A4471) 
= 0.3, EW„fe,(A6678)/EWab.(A4471) = 0.1 and EW„,ft, ( A7065) /EW„a, ( A4 471) = 0.1. The 



justification of the values of these ratios is given in llzotov fc Thuan I (120041 ) and in §4.41 The 



ionization correction factor /CF(He"'"+He^"'") is discussed in §4.6[ 



SYSTEMATIC EFFECTS 



We will consider successively the following systematic effects in the derivation of He 
abundances: 1) He I emissivities; 2) reddening; 3) the temperature structure of the H ll 
region, i.e. the temperature difference between Te(He"'") and Te(0 ill); 4) underlying stellar 
He I absorption; 5) collisional excitation of hydrogen lines; 6) the ionization structure of the 
H II region; and 7) the deviation of hydrogen and He I recombination line fluxes from case 
B. At the very end, we will also consider possible biases introduced by different samples 
not observed and reduced in the same way by deriving Yp independently for the HeBCD 
and SDSS samples and comparing the results. However, we will not discuss here collisional 
and fluo rescent enhanceraents o f the He I emission lines. These were discussed in detail 
by e.g. [Benjamin et al. I (jl999l . |2002| ) . We only note that the importance of these two 



effects depends on the specific He I emission line considered. In particular, the correction for 



-14- 



collisional excitation of the He I A5876 emission line can decrease its flux by as much as ~17% 
depending on Te(He"^) and A'"e(He"'"). In the following, we apply the correcti ons for coUisional 



and fluoresc ent enhancements, adopting the most recent analytical fits of iBenjamin et al 



(119991, [2002) 



Effects 1-6 have been considered in many previous papers on the determination of the 
primordial He abundance. However, previous analyses were typically based on small samples 
and/or not all systematic effects were considered. By contrast, the present study has two 
major advantages: 1) it is based on the largest sample of H ll regions ever assembled for the 
determination of Yp] and 2) all important known systematic effects are taken into account 
using a Monte-Carlo approach. 



4.1. He I emissivities 



First, we consider the difference in the primordial He mass fraction Yp caused by using 
two different sets of He I line emissivities. We adopt the basic set of parameters from §3.41 
except for: a) the range of TpfH e"*") variations, wh i ch we adopt here to be (0.9 - 1.0)xTe(O 
III), the same as that used by llzotov &: Thuan I (120041 ). and b) the oxygen and nitrogen 
abundances which are calculated with the electron tempera t ure s et to Te(He"'"). There are 
two differences between the procedures of llzotov &: Thuan I (120041 ) and the ones used here: 
1) the average is calculated using only three He I e r nission lines while the computation 
of takes into account all five li nes, llzotov fc Thuan I (1200411 included all five lines for the 
calculation of both y+ and 2) llzotov fc Thuan I J2004h ~adopted EWafe,(H7 + A3889) = 
3. OA instead of the relation EWafts(A3889) = EWa6s(A4471) adopted here. However, we will 
show later in this section that the variations of the EWafes(A3889) have little effect on the 
derived Yp. 

In Fig. [2^ and [2)3 we show the linear regressions Y - 0/H and Y - N/H fo r the 



HeBCD sample where the values of Y are calculated with the IBenjamin et al. I (|2002| ) He I 
emissivities. From these regressions we derive Yp = 0.2440 ± 0.0013 and 0.2464 ± 0.0010. If 
EWa6,(H7 + A3889) is set to 3.0A instead of adopting the relation EW(A3889) = EW(A4471), 
then Yp = 0.2435 ± 0.0013 and 0.2462 ± 0.0010 respectively for the y - 0/H and F - N/H 
regressions. These values are i n agre ement with Yp = 0.2421 ± 0.0020 and 0.2446 ± 0.0016 



obtained by llzotov fc Thuan I (120041 ) for their sample of 7 H ll regions. Note that the value 
of Yp derived from the Y - N/H regression is always slightly greater than the one derived 
from the Y - 0/H regression. This is b ecause the N/0 abundance ratio tends to increase 



with increasing oxygen abundance (e.g. llzotov &: Thuan Ill999l : llzotov et al. 1120061 ). 
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Figs. ^ and[2ji show the hnear regressions Y - 0/H and Y - N/H when the lPorter et al. 



(120051 ) He I emissivities are used. We obtain Yp = 0.2482 ± 0.0012 and 0.2507 ± 0.0009 from 
these regressions. It is seen that the use of the new emissivities increases 1^ by ~ 1.7%. 

In the analysi s of the other s y stema tic effects, we will consider mainly helium abundances 
obtained with the iPorter et al. I (120051) emissivities, w ith only occasional mention of helium 
abundances obtained with the iBenjamin et al. I (120021 ) emissivities for comparison. 



4.2. Reddening 



Izotov. T 



the Whitford 



luan fc Lipovetskv I (Il994l . ll997f ) and llzotov fc Thuan I (Il998al . l2004h have used 
(119581 ) reddening curve with Ry = Ay /E{B — V) = 3.2 to correct the line 



intensities for extinction. The derived value of Yp should not be very sensitive to the choice 
of a particular reddening curve because of the way the extinction correction is applied: 
the fluxes of the Balmer lines are corrected for extinction so that the resulting flux ratios 
correspond to the theoretical values, independently of the adopted reddening law. Then, the 
corrected fluxes of the other emission lines are not sensitive to the particular reddening law 
as well. 

To check the s e nsitiv ity of Yp on the adopted reddening law, we have also considered 
the ICardelli et al. I (119891) reddening curve s with Ry = 2.5, 3.2 and 4.0. We flnd that the 

(119891) redden i ng cu rve and Ry = 3.2 is only ~ 0.3% 



Yp derived with the ICardelli et al. 



higher than the value derived with the IWhitford I (11958! ) reddening curve. The Yp derived 
with Ry = 2.5 is ~ 0.4% lower than the one with Ry = 3.2, while the Yp derived with Ry = 
4.0 is ~ 0.2% higher. The effect of various reddening laws on the derived Yp is indeed small . 
Therefore, for consistency with our previous work, we have adopted the Whitford ( 1958h 
reddening la w w ith Ry = 3.2, the same as the one used by llzotov. Thuan &: Lipovetsky 



(Il994j . Il997f ) and llzotov fc Thuan I (Il998al . 120041 ) 



4.3. Temperature structure 



To derive the He abundances, llzotov. Thuan fc Lipovetsky I (119941 . 119971 ) and llzotov fc Thuan 



(Il998bl ) have assumed that the temperatures Te(He"'") and Te(0 ill), averaged over the whole 
H II region, are equal. Te(0 ill) is determined from the observed [O lll]A4363/(A4959 + 
A5007) emission line flux ratio. However, because of the high sensitivity of the flux of the 
auroral [O lll]A4363 emission line to temperature, Te(0 ill) tends to be characteristic of the 
hotter zones in the H ll region, and may be expected to be higher than Te(He"'"). On the 
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other hand, Te(H"'") and Te(He+) are derived from the recombination spectrum of ionized 
hydrogen and hehum, which is not very sensitive to temperature. Therefore, the assump- 
tion Te(He+) = Te(H"'"), which we adopt here, seems reasonable, the and He^ zones in 
high-excitation H ll regions being nearly coincident. 



To account for the difference between Te(0 ill) and T^fHe"*"). iPeimbert I (119671 ) has de- 
veloped a formalism based on the average temperature Tq and the mean square temperature 
variation in an H ll region. Then Te(He"'") and Te(0 ill) can be expresse d as functions of 
Tn an d t^, with Te(0 III) > Te(He+). This approach has been applied by iPeimbert et al. 



( 120021 ) for the determination of the He abundance in several low-metallicity dwarf galaxies. 
They find that the difference between Te(0 ill) and Te(He+) results in a reduction of the He 
mass fraction by 2 - 3 percent as compared to the case with Te(0 ill) = Te(He+). 

Until very recently, no direct measurement of Te(H+) has been carried out for extremely 
metal-deficient H ll regions, those which carry the most weight in the determination of Yp. All 
existing measurements of Te(H+) based on the Balmer jump or the Paschen jump, have been 
done for relatively metal-rich H ll regions. The lowest-metallicity H ll region for which such 
a me asurement was performed is Mrk 71 with 12 + log 0/H ~ 7.9 (j Gonzalez- Delgado et al. 
19941 ). ~ 5 - 6 times more metal-rich as compared to the lowest-metallicity H ll regions in our 
sample. Other measurements of TpfH"*") have been done for 30 Dor in the Large Magellanic 



Cloud (12 + lo g 0/H = 8.3. IPeimbert 



in ou r Galaxy (IPeimbert et al. 
2005h . 



1993 



2000 



20031) and some nearly solar-metallicity H ll regions 



Esteban et al. Ill998t iGarcia-Rojas et al. 112004 



Using the Balmer and Paschen jumps, iGuseva et al. I (120061 ) have measured Te(H"'") 
for an extensive sample of 47 H ll regions which includes the most metal-deficient BCDs 
known, such as SBS 03 35-0 52W, I Zw 18 and SBS 0335-052E. Contrary to the sugges- 
tions of IPeimbert I (119671 ) and IPeimbert et al. I (120021 ). they found no statistically significant 
difference between Te(0 ill) and Te(H"'"), but could not exclude small differences of < 5%. 

To investigate the effect of these possible small differences between Te(0 ill) and Te(He"^) 
= Te(H+) on Yp, we have derived Y - 0/H regression fits for the two sets of He I emissivities, 
both for the case where Te(He+) is varied in the range (0.95 - 1.0) x Te(0 ill) (Fig. [3^), and 
for the one where Te(He"'") = Te(0 ill) (Fig. [3)i). AH other parameters are those from the 
basic set. Comparison of Fig. [3^ (the one corresponding to the basic set of parameters) and 
Fig. [3)i shows that assuming Te(He^) = Te{0 ill) instead of varying it in the range (0.95 - 
1.0) X Te(0 III) gives a Yp larger by about ~ 1.0%. We note that if Te(He^) is different from 
Te(0 III), it may also be used for the determination of 0/H (Fig. [3}). If we do so, then Yp is 
increased by only < 0.3% compared to the basic regression in Fig. [3^, but the slope of the 
regression line is significantly shallower as compared to the case where oxygen abundance 
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is calculated with Te(0 ill). This is because the oxygen abundance derived with Te(He"'") is 
larger than the one derived with T,{0 ill), and A(0/H) = 0/H[re(He+)] - 0/H[Te(0 ill)] 
is larger for high-metallicity H ll regions than for low-metallicity ones. 

We have also co nsidered the changes o n Yp if we adopt the temperature variations of 



Te(He"'") proposed by lPeimbert et al. I (120021 ). Their Fig. 1 shows that Te(He"'") varies within 
a range (0.97 - 1) xTe(0 ill) for H ll regions with Te(0 ill) = 20000K, with the mean square 
temperature variation ranging between and 0.04, but increases to (0.88 - l)xTe(0 ill) 
for H II regions with Te(0 ill) = lOOOOK, for the same range of t^. By adopting the above 
temperature dependence of Te(He"^), we obtain Yp = 0.2506, essentially the same value as 
the one obtained from the basic regression. 



4.4. Underlying stellar He I absorption 



19861 : IPagel et al. I Il992 



Olofsson 



It has long been recognized (IRavo et al. Ill982l : iKunth fc Sargent Ill983l : iDinerstein fc Shields 



I995I ) that absorption caused by hot stars in the He I 
lines can induce an underestimation of the intensities of the nebular He I lines. In par- 



ticular, llzotov fc Thuan I (jl998al ) have shown that the neglect of He I underlying stellar 
absorption has led to the derivation of a very low helium r aass fraction in I Zw 18, the 



second most met a l-defic i ent BCD known (Pagel et al. Ill992l: Olive et al. 19971 ) . Rec ently. 



Izotov fc Thuan I (120041 1 . lOhve fc Skillman I (120041 1 and iFukugita fc Kawasaki I (120061 1 have 



taken into account thi s effect to derive Yp, using a subsample of the HeBCD sample of 



Izotov fc Thuan I mo4 ]. 



Gonzalez Delgado et al. I (119991 ) have calculated synthetic spectra of H Balmer and He 



I absorption lines in starburst and poststarburst galaxies. They predict the equivalent width 
of the He I A4471 absorption line to be in the range ~ 0.4 - 0.5 A, or < 10% of the He I 
A4471 emission line equivalent width for young starbursts with an age < 5 Myr, which is the 
case for the H ll regions in our sample. 

Unfortunately, those authors did not calculate absorption line equivalent widths for the 
other prominent He I lines. We expect however that the effect of underlying absorption is 
smaller for the He I A5876, A6678 and A7065 emission lines. This, for the following reason. It 
is known that the equivalent widths are not the same for all hydrogen absorption lines. For a 
fixed age of the stellar population, EWabs is the largest for H(5, then decreases progressively 
for H7, H/3 and Ha at longer wavelengths. A similar trend is likely to hold for He I absorption 
lines: the longer the wavelength of the line, the smaller its EWabs. It is not clear, however, 
how strong is underlying absorption for the He I A3889 line, at the short wavelength end. In 
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any case, underlying stellar absorption must be taken into account for all He I lines if we are 
to achieve the desired high accuracy of <1% in the primordial He abundance determination. 

Among the five He I emission lines used here, the effect of the underlying stellar ab- 
sorption is most important for the blue He I A4471 emission line because it has the lowest 
equivalent width. For the redder lines, especially for the He I A5876 emission line which 
carries the highest weight in the Y determination, the effect of underlying stellar absorp- 
tion is likely considerably lower. As for the blue He I A3889 emission line, its flux is > 2 
times greater that of the He I A4471 emission line and hence the effect of the underlying 
stellar absorption is also expected to be lower. However, the situation with this line is more 
complicated compared to other lines because of its blending with the hydrogen emission line 
H7. 

We have presented our choice for the various ratios linking the absorption equivalent 
width of the He I A4471 line to those of the other He I lines in our discussion of the basic 
set of parameters. While the choice of those ratios may appear arbitrary and is somewhat 
subjective since, with the exception of the A4471 line, we do not have any guidance from 
observations or models concerning the equivalent widths of the He I absorption lines, we can 
estimate the effect of different adopted EWabs on Yp. Thus, assuming larger EWabs for the 
He I A5876, A6678 and A7065 lines, in the range 0.5 - 1.0 that of the A4471 line, would result 
in a value of Yp that can be as much as ~ 1% larger than the basic value. However, our 
calculations show that for the majority of the H ll regions, the Xmm values obtained with 
the larger EWabs for the A5876, A6678 and A7065 lines are ~ 2 times larger than the Xmin 
value obtained with the basic parameter set. Additionally, the dispersion of the points in 
the Y - 0/H and Y - N/H diagrams is higher. We also made calculations with the case 
EW„6,(A5876) = 0.1xEWafc,(A4471). The resulting Yp is only <0.3% lower than the basic 
one. In summary, variations of EWafe^s for He I lines other than A4471 do not affect much 
Yp. Therefore, we will adopt the basic EWabs ratios in our Yp determination. 

Once the EWabs ratios chosen, we need to fix the absorption equivalent width of the 
A4471 line. To check the influence of the correction of He I emission line fluxes for underlying 
stellar absorption on Yp, we consider two cases: case a) where there is no underlying absorp- 
tion, i.e. EWabs(A4471) =0 and case b) where EWafes(A4471) = O.SA, the maximal equivalent 
width expected for underlying absorption. In Fig. [3)d and[3t, we show the linear regressions 
respectively for cases (a) and (b). It is seen that the difference in Yp between the two cases is 
~ 3%. In addition, we have also checked how variations of EWaf,s(A3889) affect Yp. Since the 
He I A3889 emission line is blended with the hydrogen H7 emission line, we have considered, 
in addition to the basic case shown in Fig. [3^ where EWaf,s(A3889) = EWafes(A4471), also 
the case where EWabs (H7 + A3889) = 3.0A, shown in Fig. [3li. The second case gives a Yp 
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which is ~1.0% lower than the value given by the basic model. However, it has a Xmin which 
is significantly larger than the one of the basic model. 

The effect of underlying stellar absorption is greater in H ll regions with lower equivalent 
widths of He I emission lines. This allows us to check whether the value of 0.4A which we 
have adopted for EWa{,s(A4471) in the basic set is reasonable or not. If the adopted value 
is not correct, then the values of Yp derived from H ll regions with lower EW(H/3)s and 
higher EW(H/3)s will be different, as the first ones will give underestimated He abundances 
as compared to the latter ones. In Figs. [3fe and[3]F, we show respectively the linear regressions 
for the subsamples of H ll regions with EW(H/5) > lOOA and of those with EW(H/5) > 200A. 
It is seen that the Yp derived for the HeBCD subsample with EW(H/5) > lOOA (Fig. ^) is 
only slightly (< 0.2%) lower than the value derived for the HeBCD subsample with EW(H/?) 
> 200A (Fig. [3f). This shows that EWabs(Re I A4471) = 0.4A and the ratios of EW^fes 
adopted in the basic model characterize well the underlying He I stellar absorption. 



Probably, as proposed by iPorter et al. I (120061 ) . the best check of the appropriateness 



of the basic set of EWab^s comes from comparing the abundance of singly ionized He y~^, 
calculated separately for each of the three He I A4471, A5876 and A6678 emission lines, with 
their weighted mean y^^. We show in Fig. the dependence of {y^ - yt,m) the He I 
emission line equivalent width EWem, for each of the three emission lines, in the case where 
underlying absorption is not taken into account. It is seen that there are no significant offsets 
between the y^s derived from the He I A5876 and A6678 emission lines (labeled respectively 
by stars and open circles), suggesting that the effect of underlying absorption for those 
lines is relatively small. On the other hand, the y^s derived from the He I A4471 emission 
line (filled circles) are systematically lower than those derived from the other two emission 
lines, suggesting that underlying stellar absorption is significant. Figs. - IDi show the 
same dependences as in Fig. HIi, but for successively increasing values of EWa6s(A4471), 
varying from 0.2A to O.sA. For the other He I lines we adopt the same ratios of EWah^ to 
EWab<j(A4471) as those in the basic model. It seen from Figs. Ht that the y^s derived for He 
I A4471 with EWabs(A4471) = 0.4A are in good agreement with the ones for the other two 
lines. However, the ?/"'"(A4471) derived with lower and higher EWab<j(A4471)s are respectively 
lower and higher than t he y^ values d erived for the He I A5876 and A6678 lines (compare 



Fig. Hb with Fig. Hli). iPorter et al. I (120061 ) have also suggested to use a simple mean ?/+ 
instead of the weighted mean y^^ used in our analysis. In Figs. [5^. and [Sb we show the 
dependence of (?/+ - y+) on EWem for two choices of EWafes(A4471): (Fig. [5^) and 0.4A 
(Fig. [5b). Similarly to Fig. H^, there is a clear difference between the y^ derived from the 
He I A4471 line and those derived from the other two He I lines when underlying stellar 
absorption is not taken into account (Fig. [5^). On the other hand, there is no clear offset 
in the y^ derived from the different He I lines in Fig. [5b, where EWabs(A4471) = 0.4A. The 
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primordial He abundance Yp derived from the Y O /H regression, using a simple mean for 
and EWa6s(A4471) = (Fig. [St) is ~0.8% lower than the Yp obtained using a weighted 
mean for (Fig. [Hb)- This is because the He I A4471 line contributes more to the simple 
mean than to the weighted mean. On the other hand, no difference in Yp is found when 
EWa6s(A4471) is set to 0.4A, irrespective of whether is computed with a weighted mean 
or a simple mean (compare Figs. [Sji and[3K). 

All the above considerations lead us to conclude that the basic set of EWahsS adopted 
for the He I lines is the most appropriate one to use. 



4.5. Collisional excitation of the H lines 



It has generally been assumed in abundance studies that deviations of the observed 
Ha/H/3 flu x ratio from the theor e tical recombination value are entirely due to interstellar 



extinction. 



Davidson fc Kinman I (119851 ) first noted that in the hot and dense H ll regions 



of BCDs, collisional excitation of hydrog e n em ission lines can be important and affect the 



derived He/H ratio. IStasihska fc Izotov I (120011 ) estimated that this effect can result in an 



upward correction in the He abundances of up to 5%, assuming that the excess of the Ha/H/3 
flux ratio above the theoretical recombination value is due only to collisional excitation. 



(Anderson et al. 


2002) dve 


excitation. 


Luridiana et al. 



Luridiana et al. I (120031 ) have performed a very detailed analysis of 3 H ll regions, 
including a tailored photoionization analysis and a discussion of reddening. They find that 
the collisional contribution to Ha in those objects may reach 8 % and that to H/5 2-2.5 
%. Thus, in our Monte Carlo approach, we have varied the fraction of the Ha emission line 
flux due to collisional excitation, A/(Ha)/J(Ha), in the range - 5%. Its value is derived 
together with other parameters by the minimization of &s defined in Eq. [TTl To check 
how the inclusion of hydrogen collisional excitation influences Yp, we have also calculated 
the linear regression for the case where that effect is not taken into account, i.e. AJ(Ha) = 
0. The result is shown in Fig. [3jc- Comparison with the basic regression (Fig. [3k) shows 
that ignoring hydrogen collisional excitation results in a decrease of 1^ by ~ 1.2%o 



•^In the above analysis, wc have adopted a ratio of Ha to H/3 emissivities due to colUsional excitation 
equal to 3. The referee has pointed out that this ratio may be high er. We have estimated it to be 5 - 7 in 
the temperature range Te — lOOOOK - 20000K, consistent with the lAnderson et al. I (|2002l ) determination. 
The consequence of adopting a ratio of Ha to H/3 emissivities due to collisional excitation equal to 5 would 
be to increase Yp by only 0.3%, using the HeBCD sample and in the case of the basic model. 
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4.6. Ionization structure 



Another source of systematic uncertainty comes from tlie assumption tliat tlie H"*" and 
He"*" zones in the H ll region are spatially coincident. However, depending on the hardness 
of the ionizing radiation, the radius of the He"*" zone can be smaller than the radius of the 
H"*" zone in the case of soft ionizing radiation, or larger in the case of hard radiation. In 
the former correction for unseen neutral helium should be made, resulting in an 

ionization correction factor /CF(He"'" + He^"*") > 1 and hence a higher helium abundance. 



correction factor problem has been discussed in 


several studies ( 


Stasihska 1980: Paeel et al. 


1992 


Steieman et 


al. 


1997 




Olive et al. 


1997 




Vieeas et al. 


2000: Peimbert et al. 2000: 


Ballantvne et al. 


2000 


Sauer & Jedamzik 


20021). It was concluded that the correction of 



the helium abundance can be as high as several percent in either downward or upward 
directions, dependi n g on the hardness of the radiation and the ionization parameter U. 



Sauer fc Jedamzik I (120021 ) have calculated an extensive grid of photoionized H ll region 
models which give correction factors as functions of hardness and U. Their conclusion was 
that a downward correction of Y by as much as 6% and 2% is required respectively for 
ionization parameters log U = -3.0 and -2.5. However, if log U > -2.0, which is the case for 
the majority of our H ll regions, the downward correction is < 1%. 



Ballantvne et al. I fl2000f ) have suggested that ICF{}ie+ + He2+) can be estimated from 
the [O lll]A5007/[O l]A6300 emission line flux ratio. They have calculated an extensive grid of 
photoionized H ll region models, spanning a wide r ange of metallicity and ex citation and with 
different models for the ionizing stellar radiation. iBallantyne et al. I (120001 ) have concluded 
that the ICF can be significantly lower than unity in hot H ll regions. However, if a H ll 
region has a [O ill] A5007/[O l]A6300 ratio greater than 300, then their models show that, 
regardless of its metallicit y , it will have an ICF very close to unity. Using the observations of 
Izotov fc Thuan I Jl998bh . IBallantyne et al. I (120001 ) have found that the scatter of the data 



points for H II regions with [O ill] A5007/[O l]A6300 > 300 in the Y - 0/H diagram is some 
20% smaller than that for the whole sample, and have derived Yp = 0.2489 ± 0.0030 for the 
restricted sample. However, their analysis suffers from the small number of H ll regions in 
the restricted sample. In particular, they found a negative dy/d(0/H) slope, which leads 
to a high value of Yp. Such a negative slope is unphysical and cannot be justified by any 
chemical evolution model, with any star formation scenario. 

To check how the ICF is related to the [O ill] A5007/[O l]A6300 ratio, we show in Fig. 
[6]y - O/H regression lines for the combined HeBCD + SDSS sample, with different cutoffs 
of the [O III] A5007/[O l]A6300 ratio. The HeBCD H ll regions are shown by filled circles 
and the SDSS H ll regions by open circles. We adopt JCF(He+ + He^+) = 1 for all H ll 
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regions. All other parameters are those of the basic set ( §3.41) . It is seen that the scatter in 
the Y - 0/H diagram of the HeBCD H ll regions with [O ill] A5007/[O l]A6300 > 300 (Fig. 
Eti ) is ~ 20% lower than th at for the total HeBCD sample (Fig. [B^), confirming the finding 



by iBallantyne et al. I (120001 ). However, for the SDSS sample which is larger and hence has 
better statistics, we do not find that the scatter decreases with increasing value of the cutoff. 
This suggests that the scatter of points of ~ 3% - 4% around the regression hne in Fig. 
|6] is due mainly to observat ional statistical uri c ertain ties, and not so much to variations of 
the ICF. By contrast with IBallantyne et al. I (120001 ). we find only small changes in 1^, by 
< 0.5%, when the cutoff is increased. These very small Yp changes can be due in part to 
the progressively smaller size of the H ll region sample with increasing cutoff, and hence to 
larger uncertainties in the regression. The nearly constant value of Yp implies that /C-F(He"'' 
+ He^"*") is close to the adopted value of unity for our H ll regions, regardless of the value of 
the [O III] A5007/[O l] A6300 ratio which varies in the range ~ 30 - 700. 



To derive /CF(He"'" + He^+), we use the photoionized H ll region r nodels bvlStasihs. 



(120031 ) , but with an input radiation field computed with Starbur s t 99 ( Leitherer et al. 



ca fc Izotov 



19991 ) 



using the stellar model atmospheres described in ISmith et al. I (120021 ). In Fig. [71 we show 
JCF(He"'" + He^+) as a function of the 0^+/(0"'"+0^+) abundance ratio, the latter being 
a measure of the H ll region excitation. Filled circles represent low-met allicity H ll regions 
with 12 + log 0/H = 7.2, open circles intermediate-metallicity H ll regions with 12 + log 
0/H = 7.6 and stars high-metallicity H ll regions with 12 + log 0/H = 8.3. It is seen 
from Fig. [7] that the ICF in low-excitation H ll regions, those with low 0^+/(0++0^+) 
abundance ratio, can significantly deviate from unity. However, our sample consists only 
of high-excitation H ll regions with 0^^/(0^+0^^) > 0.5 and for these objects, ICF de- 
viates from unity by not more than 1%. In our calculations, we use the ICFs derived by 
interpolating or extrapolating in metallicity the values given by the three curves in Fig. [TJ 
for a given x(0^+) = A^(0^+)/iV(0+-|-0^"*'). We find that the effect of ionization structure 
on the derived Yp is small. Neglecting the correction for ionization, i.e. adopting ICF{B.e'^ 
+ He^"*") = 1, we obtain a Yp which is only <0.5% larger than the value derived when the 
correction for ionization is taken into account (compare Fig. [31 and[3K). 



4.7. Deviations from case B recombination 

All previous work on the primordial He abundance determination has assumed that case 
B recombination holds in H ll regions. Case B assumes that there exists a balance between 
the absorption and emission of photons in the resonant Lyman series transitions of hydrogen 
and helium and that there are no other processes. However in real H ll regions, processes 
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do exist that may disrupt this balance. The first of these processes is the leakage of Lyman 
line photons from the H II region because of a finite optical depth. The second process is the 
absorption of UV photons by dust inside the H ll region. The third process is line pumping 
of stellar photons. 

Thanks to our large data base, we are able to estimate the effect of the deviations from 
case B for He I emission lines. This is because some singlet emission lines, such as the He 
I A5015 and A7281 lines, are sensitive to this effect, while other singlet lines, such as the 



He I A4922 and A6678 lines, are much less sensitive (lAller Ill984l ). Therefore, singlet He I 
emission line flux ratios may be good indicators of deviations from case B. The SDSS H ll 
region sample is particularly important in this regard, because the relatively high spectral 
resolution of the SDSS spectra allows to separate the He I A5015 emission line from the 
strong [O III] A5007 emission line. 

In Fig. El we show the dependence of the He I A5015/A6678, A7281/A6678 and A4922/A667^ 
singlet flux ratios on oxygen abundance 12 + log 0/H for the SDSS H ll regions. We divide 
the SDSS objects into two groups: the filled circles represent H ll regions with EW(H/?) > 
lOOA and the open circles those with EW(H/?) < lOOA. This division allows us to check 
whether or not deviations of the He I flux ratios from case B are affected by underlying 
stellar absorption. This effect is expected to be higher in H ll regions with low EW(H/5). 
If it plays a role, then a systematic shift between H ll regions with high and low EW(H/3) 
should be present in Fig. [HI The solid and dashed lines in Fig. [S] show respectively the 
theoretical ratios for cases B and A, in the case of two electron temperatures, = 10000 K 
and Te = 20000 K. In panels a) and c), the two lines for case A are nearly indistinguishable. 
The emission line which is most sensitive to deviations from case B is He I A5015. It is seen 
from Fig. [8^ that the He I A5015/A6678 flux ratio is systematically lower than the theo- 
retical value for case B. The same is true, but to a lesser extent, for the He I A7281/A6678 
flux ratio. It is also seen that there is no systematic offset between H ll regions with high 
and low EW(H/3). The situation is different for the He I A4922/A6678 flux ratio. There 
is no statistically significant deviation for H ll regions with high EW(H/5), but the ratios 
for H II regions with low EW(H/5) are systematically below the case B theoretical value 
(Fig. [Ht). Thus, underlying absorption does play a role, and its effect is most important 
for the He I A4922 emiss i on lin e. This conclusion is supported by the synthetic spectra of 



Gonzalez Delgado et al. I (119991 ) who predict rather high (0.3 - 0.8 A) equivalent widths for 



the He I A4922 absorption line. 

Thus, the He I A5015/A6678 and A7281/A6678 emission line flux ratios clearly show 
deviations from case B. It is reasonable to suppose that deviations from case B are also 
present for hydrogen. However, the overall effect of these H deviations on the derived element 
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abundances is likely to be small. I Cot a fc Ferland I (119881 ) have considered the effect of case 
B deviations on hydrogen line emissivities caused by a finite optical depth of the Lyman 
series lines, the presence of velocity gradients and dust. Applying their results to the low- 
metallicity H ll regions in our sample, we find that effect to be small. This, because of 
several reasons. First, the Lya optical depth in our H ll regions is likely to be very large, 
more than ~ 10^, as ev idenced by the very high H I column densit ies > 10^^ cm~^ measured 



in some of our objects (IThuan fc Izotov 1119971 : iKunth et al. Ill998l ). Second, the ratio of the 
number of dust particles to the number of hydrogen atoms in our objects is more than one 
or der of magnitude low er than the ratio in the interstellar medium of the Galaxy, as shown 
by llzotov et al. I (120061 ) in their study of SDSS extragalactic H ll reg ions. For low-metallicity 
H II regions with these characteristics, we find from Figs. 3 and 4 of I Cot a fc Ferland I (119881 ) 
that the recombination rate for hydrogen is smaller than the rate for case B by not more 
than 1%. Similar conclusions can be reached for helium. 

As for case C which o ccurs when ab sorption of stellar photons at the wavelengths of 
Lyman lines is important (IFerland I Il999l ). our Cloudy (version COG. 02) calculations show 
that this situation does not happen in any significant degree in H ll regions such as the ones 
in our sample. 

For the objects considered here, the effects of deviations from case B work in the same 
direction for both hydrogen and helium. Since the determination of the He abundance 
involves the ratio of A^(He) to A^(H), these effects partially compensate each other. In 
the present paper, we have chosen not to take into account the deviations from case B in 
the calculation of He abundances. Accounting for them accurately would require a precise 
photoionization modeling for each individual object which is not feasible, given the large size 
of our sample. However, it would be interesting to attempt this modeling for a few cases in 
the future, for checking purposes. 



4.8. Effect of different samples 



In our previous papers (llzotov. Thuan fc Lipovetsky Ill994j . 119971 : llzotov fc Thuan Ill998b 



20041 ). we have only considered H ll regions from the HeBCD sample, data which we acquired 



and reduced ourselves. In this section, we consider the effect on the Yp value when the H ll 
regions from the SDSS sample, observed and reduced independently by the SDSS team, are 
added. In Fig. [3], we show the linear regression Y - 0/H obtained with the basic parameter 
set for the combined HeBCD + SDSS sample. It is seen that SDSS data (open circles) have 
a larger scatter as compared to the HeBCD data (filled circles) which is expected because 
of the lower signal-to-noise ratio of the SDSS spectra. However, there is no visible offset 
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between the two data sets. The derived Yp value for the combined sample is close to the one 
obtained for the basic case (compare Fig. [3^ and[3l). The difference in 1^ is ~ 0.5%. 

We have also derived the regression line for the SDSS subsample only with the basic 
parameter set. We obtain Y = (0.2459±0.0021) + (90.60±13.5) 0/H, as compared to 
Y = (0.2509±0.0012) + (52.48±9.0) 0/H for the HeBCD subsample only, with the same 
parameters (see Fig. [3k). The YpS for the two subsamples are consistent with each other at 
the 2a level. The SDSS subsample was not designed specifically for the Yp problem. The 
derived lower Yp and the steeper slope are due to both the lack of very low-metallicity H ll 
regions and the lower quality spectra in the SDSS sample. 



4.9. Budget of systematic errors 



Table E] gives a summary budget of the various systematic errors discussed above. After 
the important effect of coUisional enhancement of the He I lines (not listed in Table 
which can be, depending on the particular line, larger than 10%, the second most important 
systematic effect is that of underlying He I stellar absorption, as already emphasized by 
Izotov fc Thuan I (Il998al ). Its effect is < 3%. The remaining systematic effects - temperature 
fluctuations, coUisional excitation of H lines and ionization correc tion factors - ar e all < 
1%, some increasing Yp, and some decreasing it. If we adopt the [Benjamin et al. I (l2002f ) 



emissivities as in our previous work, then all systematic effects list ed in Table [5] add up to 
a net increase of < 2%. Thus, the value Yp = 0.243 obtained by llzotov fc Thuan I (120041 ) 
for the HeBCD sample without considering the systematic effects listed in Table |5] should 
be increased to ~ 0.248. This is indeed the case. Examination of the maximum likelihood 



linear regressions in Table M for the entries with [Benjamin et al. I (120021 ) emissivities and N 



93 sho w that = . 247± 0.001, whether 0/H is calculated with Te(0 ill) or Te(He+). 



Using the lPorter et al. I (120051 ) emissivities increases Yp by <1.7%, so that Yp becomes equal 
to ~ 0.251. 



5. The primordial He mass fraction Yp and the slope dY/dZ 



To account for the various systematic errors, including those caused by coUisional and 
fluorescent enhancements of He I lines and those summariz ed in Table [5l w e have used 
a Monte Carlo method as described in §3.31 and applied by llzotov &: Thuan I (120041 ) to a 
selected sample of 7 H ll regions. We need to fix some parameters and vary others within 
"reasonable" ranges. We discuss our choices below. 
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ing from the 


Beniamin et al. 


Porter et al. 


(2005 


) changes 1 



1.7%. Because the iPorter et ah I (120051 ) emissivities 
have been calculated with updated atomic data, we will adopt them. But to compare with 
our previous work, and t o assess the change in Y p caused by the change in emissivities, we 
have also derived Yp with iBenjamin et al. I (120021 ) emissivities. As for the equivalent widths 
of He I stellar absorption lines, we adopt the basic ones ( §3.4p . This set of equivalent widths 
appears to be the most reasonable for the reasons discussed in §4.41 and also because it gives 
the lowest Xmm derived Ys in individual H ll regions. For Te(He"'"), variations in the 

range (0. 95 - 1.0) xTp(0 ill) appear to be reasonable and are consistent with the measure- 
ments of ICuseva et al. I J2006I . boOTh . We vary Ne{'iie+) and r(A3889) in the ranges 10 - 
450 cm~'^ and 0-5. The ionization correction factor /C-F(He"'"+He^+) is interpolated or 
extrapolated from the curves in Fig. [71 

In Table |3] (available only in electronic form), we show for the whole HeBCD+SDSS 
sample oxygen abundances 0/H and nitrogen abundances N/H calculated for two different 
cases, one where the electron temperature is set to Te(0 ill) and one where it is set to 
Te(He^), and He mass fractions Y. Table HI (also in electronic form only) shows the values 
of the parameters with which the abundances in Table [3] are calculated. 

Linear regre s sions F - 0/H and Y - N/H for the HeBCD sample calculated with the 



Benjamin et al. I (120021 ) He I emissivities are shown in Fig. [HI in the cases where 0/H 
= 0/H[Te(0 III)] (Fig. [Si,b) and 0/H = 0/H[T,(He+ )] (Fig. [nt,d). The corresponding 
regression lines calculated with the iPorter et al. I J2005h He I emissivities are shown in Fig. 
ITOl Comparison of Figs. [9] and [10] shows that Yp does not depend sensitively on whether 
0/H is computed with an electron temperature equal to Te(0 ill) or Te(He+), while the 
slope dF/d(0/H) is quite sensitive to the adopted T^. This is because Yp is determined 
mainly by low-metallicity H ll regions for which changes in 0/H caused by different Tg are 
small compared to the whole 0/H range for the sample. On the other hand, dy/d(0/H) 
is mainly determined by high-metallicity H ll for which changes in 0/H caused by different 
temperatures are comparable to the whole 0/H range. The same is true for the Y - N/H 
linear regressions. 



From the Y - 0/H regressions for the HeBCD sample with the [Benjamin et al. I (|2002| ) 



and [Porter et al. [ (120051 ) He I emissivities and 0/H = 0/H[Te(He "•")], we obtain respectively 
Yp = 0.2472 ± 0.0012 and 0.2516 ± 0.0011 (Fig. [S;, [TOt). The parameters of the regression 
lines are also given in Table [B The values of Yp obtained from the Y - N/H linear regressions 
(Fig. [li, [Mi, Table [6[) are higher, being res pecti vely Y^ = 0.2 4 89 ± 0.0009 and Yp = 
0.2532 ± 0.0009 for the lBeniamin et al. I J2002h and IPorter et al. I J2OO5I ) He I emissivities. 
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As noted before (e.g. Ilzotov fc Thuan I |2004J . and references therein), the systematically 
higher Yp from the Y - N/H regression is due to the nonlinear dependence of the nitrogen 
abundance on oxygen abundance. This makes the Y - N/H linear regression less reliable for 
the determination of Yp. 

For compari s on, th e lin ear regression s for t he total HeBCD + SDSS sample with the 
Benjamin et al. I (120021 ) and iPorter et al. I (120051 ) He I emissivities are shown respectively in 
Figs. [TT] and [T21 The HeBCD H ll regions are represented by filled circles and the SDSS 
H II regions by open circles . Using the Y - Q/H li near regressions, we obtain Yp = 0.2457 



± 0.0010 (Fig. [TTb ) when iBeniamin et a. 



± 0.0010 (Fig. [T2b ) when jPorter et al. I (120051 ) emissivities are used. The l^s derived for 



([20021 ) emissivities are used and Yp = 0.2505 



the HeBCD+SDSS sample are consistent within the errors with the ones derived from the 
HeBCD sample alone (compare the regression parameters in Table where the total sample 
is characterized by = 364 H ll regions, and the HeBCD sample by A^= 93 H ll regions). 
As can be seen from Figures [TT1 and [T2l the scatter about the regression lines is considerably 
larger for the SDSS sample than for the HeBCD sample. Thus since the data points entering 
in the determination of the regression lines are weighted by their errors, most of the weight in 
the Yp determination still comes from the HeBCD sample, even in the total HeBCD+SDSS 
sample. Going from 93 to 364 H ll regions does not change appreciably the la error in the 
determination of Yp-. it decreases only by 0.0001. Thus, while we show the results for the 
total sample for comparison, we consider the Yp determined from the HeBCD sample, for 
which we understand well the selection effects, as being the best value. 

Using Eq. HJ we derive from the Y - 0/H linear regressions (Table [6]) the slopes dY/dZ 
= 2.38±0.45 and 2.19 ±0.39 for the HeBCD sa mple , in the case w h ere /H is calculated with 
Te (He^), and for the lBenjamin et al. I J2002h and lPorter et al. I J2005h He I emissivities re- 
spectively (Fig. [TTb . inb). These slopes are steepe r than the ones pre dicted by closed-box 



Larsen et al~ll200lh . It is reasonable to 



chemical evolution models of dwarf galaxies (~ 1, 
expect that some fraction of chemical elements, especially oxygen, escapes the galaxy due to 
e.g. supernova explosions. The slopes dY/dZ derived from the Y - N/H linear regressions 
(Fig. UmiZld) are not as s teep: dY/dZ = 1.46 ±0.25 and 1.34±0.24 respectively for the 



Benjamin et al. I (12002! ) and jPorter et al. I (120051 ) He I emissivities, and for the case where 



0/H is calculated with Te(He"'"). Again, the difference is due to the nonhnear dependence of 
the N abundance on O abundance, as discussed above. As for Yp, this makes the determi- 
nation of dY/dZ from the Y - N/H regression not as reliable as the one from the Y - N/H 
regression. The dY/dZ slopes are steeper when Te(0 ill) is used for the 0/H determination. 
In the latter case, we obtain for the HeBCD sample, from the Y - 0/H regressions, dY/dZ 
= 2.94±0.51 and 2.88±0.50 respectively for the two sets of emissivities (Fig. [TTk and [T2b) . 
The true electron temperature is probably somewhere between Te(He+) and Te(0 ill), so 
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we expect the true s lopes to be in the range 2.2 - 2.9 for both iBenjamin et aL 
Porter et aL I (120051 ) emissivities. 



fl2002h and 



Our Yp range of ~ 0.001 is much smaller than the one of ~ 0.009 claimed bvlOlive fc Skillman 



(2004). This difference is due to 1) different methods used in this paper and in lOlive fc Skillman 



(120041 ) and 2) much larger samples used in this paper (93 observations in t he HeBCD sample 



and 3 64 observations in the total sample as compared to 7 objects in the lOlive &: Skillman 



(I2OO4J ) sample). We emphasize that the use of large sample significantly reduces the uncer- 
tainties in Yp to the level that is acceptable for cosmological implications. 



6. COSMOLOGICAL IMPLICATIONS 



6.1. The baryonic mass density 



We now investigate whether our derived values of Yp are consistent with the predictions 
of SEEN and whether the baryonic mass density corresponding to Yp agrees with the one 
derived from measurements of the CME. In Fig. [121 we show by solid curves the dependence 
of the abundances of the three light el ements He, D ari d ''Li on the baryon-to-photon number 
ratio T] as predicted by SEEN (e.g. ISteigman II2OO6I ). We do not discuss ^He because the 



derivation of its primordial value from its presently observed value in the Galaxy is compli- 
cated by uncertain effects of chemical evolution. The sol id and dashed vertica l lines show the 



value of 77 and its la deviations as obtained by WMAP (ISpergel et al. II2OO6I ). The presently 



best observational abundances of the three light elements with their la deviations are shown 
by boxes. Concerning He, we show in the left and right panels the Yp derived from regres- 
sion fitting of the HeECD sample and with 0/H obtained respectively with the electron 
temperature set equal to T^^O ill) and Te(He"'"). We shall discuss the results concerning our 
preferred values of Yp, the ones corresponding to = Te(He+). With Yp = 0.2472±0.0012 
and 0.2516±0.0011 (solid boxes in Fig. [T51i) for the two different sets of He I emissivities, 
and with an equivalent number of light neutrino species equal to 3, the SEEN model gives 



and 8.7^Jo, respectively, where the error bars denote la errors. These 



r]io = Wr] = 5.5 

values cor respond to a bar y onic m ass fractions Qhli 



-0.004 



0.0201^:^^^ and 0.032l5:55| respec- 



tively for Eenjamin et al. ( 20021) an d Porter et al. (2005) He I emissivities. The value 



derived with the old lEenjamin et al. I Jiooi 

He I emissivities is much higher than Qhh'^ de- 



2006 



rived from the primor dial ^Li abundance (lEonifacio fc Molaro 1119971 : lAsplund et al. 
Eonifacio et al. 1120071 ) but is in excellent agreement with Qbh'^ = 0.021±0.002 derived from 



recent measurements of the deuterium abundance in damped Lya systems (IKirkman et al. 
2OO3I : lO'Meara et al. Il2006h and fifc/i^ = 0.0223±0.0008 derived from measurements of the 



fluctuations of the microwave background radiation by WMAP (jSpergel et al. 2006). An- 



- 29 - 



other way of say i ng th e same thing is that the observational value of Yp derived with the 
Benjamin et al. I (120021 ) He I emissivities is in excellent agreement with the theoretical SBBN 
value of 0.248. 



On the other hand, the Yp derived with the new iPorter et al. I (120051 ) He I emissivities 
is larger than the predicted SBBN value and is only consistent with it at the 2a level. 
Consequently, the corresponding = 0.0321q;qq4 is also consistent only at the 2a level 
with the value inferred from the D abundance measurements and the WMAP data. In 



principle. thelPorter et al. I (120051 ) He I emissivities should be more reliable than those from 
Benjamin et al. I J2002h . since they are based on new updated atomic data for He. 



Steigman I (120061 ) has shown that, while He may not be as sensitive a baryometer as 



D, it is an excellent chronometer (in the sense that it is sensitive to small deviations from 
the standard Hubble expansion rate) and/or leptometer (in the sense that it is sensitive to 
any asymmetry between the numbers of neutrinos and antineutrinos). We investigate next 
possible small dev iations from SBBN i f we take at face value the relatively high value of Yp 
obtained with the iPorter et al. I (120051 ) He I emissivities. 



6.2. Deviations from SBBN 



Deviations from the standard rate of Hubble expansion in the early Universe can be 
caused by an extra contribution to the total energy density (for example by additional flavors 
of active or sterile neutrinos) which can conveniently be parameterized by an equivalent 
number of neutrino flavors N^. Combining Qbh"^ = 0.00223 ± 0.0008 obtained by WMAP 



(ISpergel et al. 1120031 ) with Yp = 0.2516 ± 0.0011, we obtain ~ 3.2 (IWalker et al. Ill99lh . 



We use the stati stical x technique with the code described by iFiorentini et al. I (Il998l ) 



and lLisi et al. I (119991 ) to analyze the constraints that the measured He, D and ^Li abundances 
put on ?7 an d N,y. For the pr i mord ial D and ^Li abunda n ces, w e use respectively the values 



obtained bv lO'Meara et al. I J2006h an d lAsplund et al. I J2006h . With Yp = 0.2472± 0.0012 



Benjamin et al. II2002I ). the minimum Xmin = 9.03 6 is obtained when ?7in 



(emissivities by 

5.79 and = 2.972. With Yp = 0.2516±0.0011 (emissivities by IPorter et al. Il2005h . the 



mmimum Xn 



9.334 is obtained when rjiQ = 5.97 and N^, = 3.280. 



The joint fits of t] and N^, with the YpS derived with the two sets of emissivities are 
shown respectively in Figure [H^ and [TDo. The la (x^ - Xmm = 1-0) 2a (x^ - Xmin 
= 2.71) deviations are shown respectively by the thin and thick solid lines. We find the 
equivalent number of light neutrino species to be respectively in the range N^, = 2.97±0.16 
(2a) (Fig. [HJi), and = 3.28±0.16 (2a) (Fig. [Hd). The first value of obtained with 
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Yp derived with the ol d set of He I emiss ivities is entirely consistent with the experimental 
value of 2.993±0.011 (ICaso et al. I Il998l ) shown by the dashed line. On the other hand, 
the second value of obtained with Yp derived with the new set of He I emissivities is 
significantly higher than the experimental value, suggesting a slight deviation from SBBN. 
We note that the primordial helium abundance sets very tight constraints on the effective 
number of neutrino species. These constraints are much tighter than those derived using 
the CMS and galaxy clustering power spectra. For example, using these two sets of data, 
Ichikawa et al. I (120061 ) derive 0.8 < Ny < 7.6 at the 95% confidence level. 



Alternatively, deviations from the SBBN model can be checked by introducing two 
parameters, the ex pansion rate param eter S, and the electron neutrino asymmetry parameter 
^e, as suggested by Steigman ([2005). If deviations from the standard expansion are caused 
by an extra non-standard energy density component, then it is convenient to express that 
extra energy component in terms of an extra effective number of neutrino species AiVi, 
defined as: Ny = ?> + ANy. Then S is related to AN^, by: 



S 



1/2 



(13) 



For the SBBN model, AiV^ = 0, A^^ = 3 and 5 = 1. 



As for the parameter C,e, it is related to the difference between the number of electron 
neutrino species n,^^ and the number of electron antineutrino species rip^ by : 



^ 1.33 



n,. 



(14) 



where is the number of photons. A non-zero would result in a different n/p ratio during 
the BBN period, therefore changing the light element abundances. In particular, if > 0, 
there are more neut rinos than antine utrinos, and reactions such as n + z/g — > p + e~ drive 
down the n/p ratio (jSteigman II2005I ). 



Following ISteigman I (120051 ). we show in Fig. [T5k the S - r] diagram with the nearly 
orthogonal isoabundance curves of D (dashed lines) and He (solid lines), in the case where 
= 0. T he open circle co r respo nds to the value of the primordial D/H abundance ratio 
derived by lO'Meara et al. I (120061) and that of the primordial He mass fraction derived in 
this paper, using the iBenjamin et al. I (119991 ) He I emissivities. T he filled circle corres ponds 
to the same D/H abundance ratio, but to Yp derived with the iPorter et al. I (120051 ) He I 
emissivities. From Fig. [ISi, we obtain S = 0.996 ± 0.007 (la) and 1.020 ± 0.006 (lex) for 
Yp obtained respectively with the old and new He I emissivities. As before, we find that 
the first value is consistent within the errors with the SBBN value, 5* = 1. On the other 
hand, the second value indicates a slight deviation from the SBBN and implies the presence 
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of additional neutrino species. The same conclusion can be reached with Fig. [T5b . It shows 
the ^ V diagram with the isoabundances curves of D (dashed lines) and He (solid lines) , 
i n the case where 5 = 1. Again, it can be seen that the Yp derived with the lBenjamin et al 
He I emissiv ities is consistent with SBBN (,^e = 0), while the Yp derived with the 



Porter et al. I (120051 ) He I emissivities gives a negative value of ^e, indicating deviation from 
SBBN. 



7. SUMMARY AND CONCLUSIONS 

We present in this paper the determination of the primordial helium mass fraction Yp 
by linear regressions of a sample of 93 spectra of 86 H ll low-metallicity extragalactic regions 
(the HeBCD sample). This sample is one of the largest and most homogeneous data sets 
in existence for the determination of Yp. For comparison and to improve the statistics for 
investigating systematic effects in the determination of 1^, we have also considered a sample 
of 271 low-metallicity H ll regions selected from the Data Release 5 of the Sloan Digital Sky 
Survey (the SDSS sample). 

In the determination of Yp, we have considered several known systematic effects. We 
have used Monte Carlo methods to take into account the effects of collisional and fluorescent 
enhancements of He I recombination lines, of collisional excitation of hydrogen emission lines, 
of underlying stellar He I absorption, of the difference between the temperature Te(He^) in 
the He+ zone and the temperature Te(0 ill) derived from the [O lll]A4363/(A4959+A5007) 
flux ratio, and of the ionization correction factor ICF{lle'^ + He^"^). We have also considered 
the effects of different sets of He I line emissivities and of different reddening laws. We discuss 
the effect of possible deviations of He I and H emission line intensities from case B. 

We have obtained the following results: 

1. After the effect of collisional enhancement of the He I lines, the second most important 
systematic effect comes from underlying He I stellar absorption (~ 3%). Other effects such 
as variations in temperature, collisional excitation of hydrogen emission lines or ionization 
corrections are smaller (~ 1%). Because those systematic effects can work in either direction, 
either increasing or decreasing Yp, they tend to cancel each other, so that the net effect of 
all the con sidered systemati c effec ts is an increase of 1^ of ~ 2% as comp ared to the value 
derived bv llzotov fc Thuan lJ200i). With the old set of He I emissivities by lBenjamin et al. 



fl2002h used by llzotov fc Thuan I tooi ) , and with the electron temperature set, not to Te(0 



III), but to Te(He"'"), we derive Yp = 0.2472±0.0012 which, according to SBBN, corresponds 
to a baryonic mass fraction Qhh^ = 0.020l|]:[|[]^. This value of Q^h"^ is consistent with the ones 
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derived from deuterium abundance observations and WMAP microwave background radia- 
tion fluctuation measurements, and all three measurements concur to support the validity 
of the S BBN model (SBBN p redicts Yp = 0.248 when Q^h"^ is set to the value derived by 
WMAP, ISpergel et al. II2006I ). We have checked that the derived Yp does not depend on the 
particular sample of H ll regions used. We have also performed regression fits for the SDSS 
sample observed and reduced by the Sloan team, and found that the derived Yp is consistent 
at the 2a level with the value obtained from the HeBCD sample. 

2. On the other hand, if the new set of emissivities, based on updated He I atomic 
data, by lPorter et al. I (120051 ) is used, then we obtain Yp = 0.2516±0.0011, corresponding to 



= 0.0321q;qq4, significanly larger (at the 2a level) than the ^^h'^ values derived from 
the deuterium abundance and microwave background radiation fluctuation measurements. 
If we take the higher value of Yp at its face value, then this would imply the existence of 
small deviations from SBBN. In order to bring the high value of Yp into agreement with 
the deuterium and WMAP measurements, we would need an equivalent number of neutrino 
flavors equal to 3.303 instead of the canonical 3. 

3. The dY/dZ slopes derived from the Y - / H lin ear re gressions for t he He BCD 



sample, using He I emissivities by iBenjamin et al. I f l2002h and IPorter et al.1 J2005h 



are 



res pectively equal to 2.38 ± 0.45 and 2 .1 9 ± 0.39, consistent with previ ous determinations 
bv Izotov. Thuan &: Lipovetskv I Jl997h . Ilzotov fc Thuan I Jl998bl . EooJ using BCDs, and 
by iJimenez et al. I (120031 ) from nearby K dwarf stars. 



We have considered here, as best as we can, all known systematic uncertainties that may 
affect the determination of the primordial He abundance. However, the real situation may be 
more complicated. The most important issue that we do not yet fully understand appears to 
be the temperature structure of H ll regions. For example, detailed photoionization models 
of H II regions often underpredict the electron temperature that is measured frorn the [O 
III] A4363/(A4959 + A5007) line ratio (see e.g. references in iLuridiana et al. II2003I ). Also, 
there are indications that t he matter inside H ll regions is not chemically homogeneous 
( iTsamis fc Pequignot II2005I . Stasihska et al, in preparation), so that the oxygen abundance 
obtained by traditional H II region abundance analysis may be biased. How much this affects 
the He/H abundance ratio, as well as the slope i n the Y vs 0/H relation is not known yet. 
Finally, as suggested long ago by iBond et al. I (119831 ). very massive primordial stars may 
produce a significant amount of helium without overproducing metals, in which case the 
primordial helium abundance could be smaller than obtained from a linear extrapolation 



of th e Y - 0/H relation. Recent computations of Population HI star yields (IMarigo et al. 



20031 ) associ ated with measuremeiits of the near-infrared background suggest that this may 
be the case (ISalvaterra &: Ferrara II2003I ). 
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Fig. 1. — [O l] A6363/A6300 emission line flux ratio vs oxygen abundance. The HeBCD 
and SDSS galaxies are shown respectively by fllled and open circles. The solid line shows 
the mean value of the [O l] A6363/A6300 emission line flux ratio and the dashed lines la 
deviations from the mean. 
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Fig. 2. — Linear regressions of the helium mass fraction Y vs. oxygen and nitrogen abun- 
danc es for H ll region s in the HeB CD sample. The He I er aissivities in ( a ) and (b) are 
from iBenjamin et al. I (119991 . |2002| ) and in (c) and (d) from iPorter et al. I (120051 ). In all 
panels, Y was derived by minimizing x 

and adopting EWa6,(A4471) = 0.4 A, EWafe,(A3889) 
= EWa6.(A4471), EWafe,(A5876) = 0.3xEWafe,(A4471), EW,b.(A6678) = EWafe,(A7065) 
=0.1xEWabs(A4471). The electron temperature Te(He^) is varied in the range (0.90 - 
1.00)xTe(O III). The oxygen and nitrogen abundances are calculated setting the electron 
temperature equal to Te(0 ill). 
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Fig. 3. — Linear regressions of the helium mass fractio n Y vs. oxy g en ab undances for 
different parameter sets. The He I emissivities are from iPorter et ah I (120051 ). (a) Ys are 
derived for the HeBCD sample adopting EWabs(A4471) = 0.4 A. The electron temperature 
Te(He"'") is varied in the range (0.95 - 1.0)xTe(O ill) (basic model), (b) same as in (a), but 
for EWabs(A4471) = 0.0 A. (c) same as in (a), but for EWab5(A4471) = 0.5 A. (d) same as 
in (a), but for EWabs(H7 + A3889) = 3.0 A instead of EWab,(A3889) = EWab,(A4471). (e) 
same as in (a), but only objects with EW(H/3) > lOOA are shown, (f) same as in (a), but 
only objects with EW(H/?) > 200A are shown, (g) same as in (a), but Te(He"'") is varied 
in the range (0.9 - 1.0)xTe(O ill), (h) same as in (a), but with Te(He+) = Te(0 ill), (i) 
same as in (a), but oxygen abundance 0/H is derived adopting Te(He"'") instead of Te(0 ill), 
(j) same as in (a), but adding the H ll regions from the SDSS. (k) same as in (a), but for 
zero coUisional excitation of hydrogen lines. (1) same as in (a), but the ionization correction 
factor JCF(He+ + }ie^+) is set to 1. 
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Fig. 4. — {y~^ - y^^) derived for each of the three He I emission hnes A4471 (filled circles), 
A5876 (stars) and A6678 (open circles) vs the equivalent width EWem oi the same lines, 
is the weighted mean of the of each individual fine. Four different values of EWabs(A4471) 
have been adopted: a) EWa6,(A4471) = 0, b) EWa6,(A4471) = 0.2A, c) EW„6,(A4471) = 0.4A 
and d) EWa65(A4471) = 0.5A. Only data for the HeBCD sample are shown. 
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Fig. 5. — a) Same as in Fig. HIi, except that is the simple mean instead of the weighted 

mean ?/+^; b) Same as in Fig. Ht, except that is the simple mean instead of the weighted 

mean y^^; c) Same as in Fig. [St, except that is the simple mean instead of the weighted 

mean y^^; d) Same as in Fig. [3^, except that is the simple mean instead of the weighted 
mean y+^. 
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Fig. 6. — Linear regressions of the helium mass fraction Y vs. oxygen abundances for for 
H II regions with differ ent values of t h e [O lll]A5007/[O l]A6300 emission flux ratio. The He 
I emissivities are from iPorter et al. I (120051 ). The ionization correction factor ICF{B.e^ + 
He^"*") is set to 1. In panel a) all H ll regions are shown, while in panels b), c) and d) are 
shown H II regions with the [O lll]A5007/[O l]A6300 emission flux ratio respectively greater 
than 100, 200 and 300. 
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Fig. 7. — Ionization correction factor JCF(He+ + He^"*") as a function of the doubly ionized 
oxygen abundance fraction 0^"'"/(0"'"+0^"'") 
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Fig. 8. — Flux ratios of singlet He I emission lines vs oxygen abundance for the SDSS H 
II regions. Filled circles are H II regions with EW(H/3) > lOOA and open circles are H II 
regions with EW(H/3) < lOOA. The flux ratios for cases A and B are shown respectively by 
dashed and solid horizontal lines. 
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Fig. 9. — Linear regressions of the helium mass fraction Y vs. oxygen and nitrogen abun- 
danc es for H ll regions from the He BCD sample. The Ys are derived with the He I emissivities 
from lBenjamin et al. I (1l999l . l2002l ). The electron temperature Te(He^) is varied in the range 
(0.95 - l)xTe(0 III). The oxygen abundance is derived adopting Te(0 ill) in a) and b) and 
Te(He+) in c) and d). 
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Fig. 10. — Same as in Fig. [9] but Ys are derived with the iPorter et al. I (120051 ) He I 
emissivities. 
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Fig. 11. — Same as in Fig. [9|, but for the HeBCD + SDSS sample. HeBCD H ll regions are 
shown by filled circles and SDSS H ll regions by open circles. 
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Fig. 12.— Same as in Fig. [TOl but for the HeBCD + SDSS sample. HeBCD H ll regions 
are shown by filled circles and SDSS H ll regions by open circles. 
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a) *lie 0/H=0/H[T (0 III)] 




Fig. 13. — Dependence of the helium mass fraction (a),(d), deuterium (b),(e) and ^Li (c),(f) 
abundances on the photon-to-baryon number ratio rj. Sohd curves are predictions of the 
SBBN, sohd and dash ed vertical lines indicate the rj and la deviations derived from WMAP 
(jSpergel et al. II2006I ). Boxes show the observed light element abundances along with their 
la deviations. The solid boxes correspond to the combined HeBCD + SDSS sample and 
the dotted boxes to the HeBCD sample only. Boxes in (a) are calculated adopting 0/H = 
0/H[Te(0 III)], and boxes in (d) are calculated adopting 0/H = 0/H[re(He+)]. 
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Fig. 14. — Joint fits to the baryon-to-photon number ratio, log ?7io, and the equiva- 
len t number of hgh t neut rino species A^,^, using a analysis with the code developed 



by iFiorentini et al. 



derived wi th the 



19981) and iLisi et al. I (119991 ) (a) for the primordial abundance value 



abundance value derived with the iPorter et al. 



Benjamin et al. I (120021) He I emissivities (this paper), (D/H)p from 



O'Meara et al. I (l2006l ) and (^Li/H)p fr omlAsplund et al. I (20061) and (b) for the primordial 



(120051 ) He I emissivities (this paper). 



and the same (D/H)p and (^Li/H)p as in (a). Thin and thick solid lines ar e respectively la 
and 2cr deviations. The experimental value = 2.993 (jCaso et al. Ill998l ) is shown by the 
dashed hne. 




Fig. 15. — (a) Diagram showing the expansion rate factor S vs the baryon-to-photon number 
ratio rj. Isoabundance curves for D and He are shown respectively by dashed and sohd 
hues and labeled by t he abundance va l ues. The filled circle corresponds to the primordial 
D/H abundance from lO'Meara et al. I ( 20061) and the primordial He mass fraction derived 
in this paper with the Porter et al. f 2005h emissivit i es. T he open circle corresponds to 



the primordial D/H abundance from O'Meara et al. f 2006h and the primordial He mass 
fraction derived in this paper with the iBenjamin et al. I (119991 ) emissivities. Error bars are 
la deviations, (b) Same as in (a) except the diagram shows the electron neutrino asymmetry 
parameter vs the baryon-to-photon number ratio 7]. 
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Table 1. Emission line fluxes 



Object^- 


[O 


nl 


3727 


H12 3750 


Hll 3771 


HIO 3797 


H9 3835 


H8+He I 3889 


H6 4102 


H7 4340 




[O 


in] 


4363 


He I 4471 


He II 4686 


up 4861 


[O III] 4959 


[O III] 5007 


He I 5876 


Ha 6563 




[N 


II] 


6583 


He I 6678 


[S II] 6717 


[S II] 6731 


He I 7065 


[O II] 7320 


[O II] 7330 





CGCG 007-025 (No. 


1) 


88.7± 


1, 


.4 


2. 


.8± 


1, 


.0 


3. 


,5± 


1, 


.0 


4. 


,3± 


1 


.0 


6, 


,1± 


1. 


,0 


17. 


,6± 


1, 


,0 


23, 


.7± 


1, 


.1 


47. 


,0± 


1, 


,1 






11.9± 


1, 


.0 


3, 


.5± 


1, 


.0 


1, 


,3± 


1, 


.0 


100. 


,0± 


1 


.0 


185, 


,8± 


2, 


,1 


564, 


,2± 


5, 


,7 


12 


.2± 


1, 


.0 


349, 


,9± 


3, 


,7 






4.1± 


1, 


.0 


3, 


,6± 


1, 


.0 


9, 


,0± 


1, 


,0 


6, 


,9± 


1 


.0 


4, 


,0± 


1, 


,0 


1, 


,6± 


1, 


,0 


1, 


,3± 


1, 


,0 










CGCG 007-025 (No. 


2) 


135.2± 


2, 


.1 


2, 


,3± 


1, 


.2 


3, 


,4± 


1, 


.2 


4, 


,9± 


1 


.2 


6, 


,8± 


1, 


,2 


16, 


,7± 


1, 


,2 


24, 


.9± 


1, 


.1 


45, 


,7± 


1, 


,2 






10.2± 


1, 


.1 


3, 


,4± 


1, 


.1 


0, 


,8± 


1, 


.1 


100. 


0± 


1 


.1 


151, 


,2± 


1, 


,9 


463, 


,4± 


4, 


,8 


11, 


.1± 


1, 


.1 


325, 


,7± 


3, 


,5 






5.3± 


1, 


.1 


3, 


.3± 


1, 


.1 


13, 


,3± 


1, 


.1 


9, 


,2± 


1 


.1 


3, 


,0± 


1. 


,1 


























HS 0029-t-1748 




72.0± 


1, 


.4 


1, 


A± 


1, 


.1 


1, 


,6± 


1, 


.1 


3, 


,0± 


1 


.1 


3, 


,9± 


1, 


,1 


12, 


,7± 


1, 


,1 


18, 


.7± 


1, 


.1 


41. 


,2± 


1, 


,1 






2.5± 


1, 


.0 


3, 


,9± 


1 


.0 










100. 


0± 


1 


.0 


134, 


,1± 


1, 


,7 


413, 


,4± 


4, 


,3 


13 


,0± 


1, 


,0 


347. 


,7± 


3, 


,6 






23.7± 


1, 


.0 


4, 


,1± 


1, 


.0 


18, 


,7± 


1, 


,0 


15, 


,1± 


1 


.0 


3, 


,5d= 


1, 


,0 


3, 


,2± 


1, 


,0 


2, 


,8± 


1, 


,0 










HS 0111-1-2115 




240.9± 


3, 


.5 


































11, 


,0± 


1, 


,5 


18, 


,1± 


1, 


,5 


38, 


,6± 


1, 


,5 






7.1± 


1, 


.1 


3, 


.5± 


1, 


.1 


0, 


,8± 


1, 


.1 


100. 


Oit 


1, 


.0 


200, 


,3± 


2, 


,3 


609, 


,6± 


6, 


,2 


13, 


.0± 


1, 


.1 


381, 


,8± 


4, 


,0 






7.4± 


1, 


.1 


4, 


.1± 


1, 


.1 


13, 


,8± 


1, 


.1 


10. 


,3± 


1 


.1 


3, 


,3± 


1. 


,1 


2. 


,0± 


1, 


,1 


1, 


.4± 


1, 


.1 










HS 01224-0743 




66.7± 


1, 


.4 


2, 


.0± 


1, 


.1 


2, 


,8± 


1, 


.1 


3. 


,6± 


1 


.1 


6, 


,0± 


1. 


,1 


16. 


,6± 


1, 


,1 


23, 


.1± 


1, 


.1 


45. 


,1± 


1, 


,2 






3.5± 


1, 


.4 


3, 


.7± 


1, 


.5 










100. 


0± 


1 


.3 


153, 


,3± 


2, 


,2 


476, 


,4± 


5, 


,1 


15, 


.1± 


1, 


.4 


413. 


,5± 


4, 


,7 






27.5± 


1 


.4 


4, 


,5± 


1, 


.4 


36, 


,1± 


1, 


,5 


24, 


,7± 


1 


.5 


4, 


,4± 


1, 


,5 



























(Abridged. Table with all H II regions will be published in the online edition of ApJ) 



^Thc first 93 entries are HeBCD H II regions with names in alphabetical order. The remaining 271 entries are SDSS H II regions. Their 
names are in the format xxxxx-yyyy-zzz, where xxxxx is the middle Juhan date (MJD) of the observation, yyyy is the plate number, and zzz 
is the fiber number. 



Table 2. Emission line EWs 



Object'^ 


[O 


III 


3727 


H12 3750 


Hll 3771 


HIO 3797 


H9 3835 


H8-t-He I 3889 


H<5 4102 


H7 4340 




[O 


III] 


4363 


He I 4471 


He II 4686 


H/3 4861 


[O III] 4959 


[O III] 5007 


He I 5876 


Ha 6563 




[N 


II] 


6583 


He I 6678 


[S II] 6717 


[S II] 6731 


He I 7065 


[O II] 7320 


[O II] 7330 





CGCG 007-025 (No. 


1) 


127, 


,8+ 


0, 


.6 


4, 


,1+ 


0, 


,3 


5, 


,1+ 


0, 


,3 


6, 


,4± 


0, 


,3 


9, 


,1+ 


0, 


,2 


26 


.3+ 


0, 


,3 


42, 


,8+ 


0, 


.2 


100, 


,4+ 


0.3 






25, 


,7± 


0, 


,2 


8, 


,1± 


0, 


.2 


3, 


,4± 


0, 


,2 


270, 


,2± 


0, 


.5 


535, 


,6± 


0, 


,6 


1633, 


,0± 


1, 


1 


47, 


,4± 


0, 


,3 


1544, 


,0± 


1.3 






15, 


,7± 


0, 


,2 


15, 


,8± 


0, 


.3 


40, 


,4± 


0, 


,3 


31, 


,3+ 





.3 


20, 


,4± 


0, 


,4 


9, 


,1± 


0, 


,4 


7, 


,6± 


0, 


,4 








CGCG 007-025 (No. 


2) 


131, 


,0± 


1, 


.1 


2, 


,4+ 


0, 


.6 


3, 


,6+ 


0, 


,5 


5, 


.2+ 





.5 


7, 


,4± 


0, 


,4 


21, 


.7+ 


0, 


,8 


32, 


,8+ 





.4 


72, 


,3± 


0.5 






16, 


,5± 


0, 


.4 


5, 


,9± 





.4 


1, 


,6± 


0, 


,3 


204, 


,3± 





.8 


308, 


,8± 


0, 


,9 


952, 


,1± 


1, 


,5 


34, 


,0± 





.6 


1152, 


,0± 


2.4 






18, 


,9± 


0, 


.6 


12, 


,1± 


0, 


.6 


47. 


,2± 


0, 


,8 


34, 


,2± 


0, 


.8 


12. 


,5± 


0, 


,8 
























HS 0029+1748 




106, 


,3± 





.9 


2, 


,1± 





.6 


2, 


,3± 


0, 


,5 


4, 


.1+ 





.4 


5, 


,3± 


0, 


,3 


19 


.2+ 


0, 


,3 


24, 


,3± 





.3 


54, 


,5± 


0.3 






2, 


,9± 





.1 


4, 


,8± 


0, 


.1 










153, 


,0± 


0, 


.2 


215, 


,9± 


0, 


,2 


681, 


,9± 


0, 


,4 


30, 


,3± 





.1 


919, 


,1± 


0.5 






62, 


,4± 


0, 


,2 


11, 


,6± 


0, 


.1 


51, 


,6± 


0, 


,2 


41, 


,7± 


0, 


,2 


10, 


,6± 


0, 


,2 


10, 


,2± 


0, 


,2 


8, 


,7± 


0, 


,2 








HS 0111+2115 




117, 


,8± 


1, 


,0 


































5, 


,1± 


0, 


,4 


8, 


,8± 


0, 


,4 


21, 


,4± 


0.4 






9, 


,3± 


0, 


.2 


4, 


,7± 





.2 


1, 


,0+ 


0, 


,2 


149, 


,4± 


0, 


.4 


298, 


,6± 


0, 


,6 


912, 


,1± 


0, 


,9 


22, 


,8± 


0, 


.3 


793, 


,9+ 


1.1 






15, 


,5± 





.3 


8, 


,7± 


0, 


.3 


29. 


,5± 


0, 


,4 


22, 


,0± 


0, 


.3 


7, 


,8± 


0, 


,4 


5, 


.1+ 


0, 


,5 


3, 


,4± 


0, 


.5 








HS 0122+0743 




103, 


,1± 





.8 


3, 


,1± 


0, 


.4 


4, 


,3± 


0, 


,4 


5, 


.5+ 





.4 


9, 


,1± 


0, 


,3 


24, 


,2± 


0, 


,4 


38, 


,7± 


0, 


,4 


83, 


,6± 


0.5 






1, 


,8± 


0, 


.3 


2, 


,0± 


0, 


,3 










63, 


,0± 


0, 


,5 


98, 


,3± 


0, 


,6 


307, 


,0± 


0, 


,8 


11, 


,5± 


0, 


.4 


357, 


,9± 


1.5 






23, 


,9± 


0, 


,6 


3, 


,9± 


0, 


,5 


32, 


,4± 


0, 


,7 


22, 


,1± 


0, 


,6 


4, 


,2+ 


0, 


,6 

























(Abridged. Table with all H II regions will be published in the online edition of ApJ) 



^The first 93 entries are HeBCD H II regions with names in alphabetical order. The remaining 271 entries are SDSS H II regions. Their 
names are in the format xxxxx-yyyy-zzz, where xxxxx is the middle Julian date (MJD) of the observation, yyyy is the plate number, and zzz 
is the fiber number. 



Table 3. Oxygen, nit roKen and hel i um ab undances for the best solutions with the 

Porter et al. I (120051 ) He I emissivities 



Object^- 






O/H'^.d 


N/H'''^ 


N/H<=''= 


Y 


EW(H/3) 


CGCG 007-025 (No. 


1) 


6.0± 0.1 


6.0+ 0.1 


13.5+ 0.3 


14.0+ 0.3 


0.2493+ 0.0060 


270 


CGCG 007-025 (No. 


2) 


5.5± 0.2 


6.1+ 0.2 


12.6+ 0.5 


13.8+ 0.5 


0.2511+ 0.0069 


204 


HS 0029+1748 




11. 3± 0.4 


13.3+ 0.4 


40.0+ 1.4 


43.5+ 1.4 


0.2525+ 0.0043 


149 


HS 0111+2115 




20.7± 2.8 


20.9+ 2.8 


81.0+ 8.7 


85.0+ 8.7 


0.2697+ 0.0100 


63 


HS 0122+0743 




4.0± 0.1 


4.4+ 0.1 


9.8+ 0.4 


11.1+ 0.4 


0.2536+ 0.0055 


232 


(Abridged. Table with all H II regions will be published in the online edition of ApJ) 



^The first 93 entries are HeBCD H II regions with names in alphabetical order. The remaining 271 entries 
are SDSS H II regions. Their names are in the format xxxxx-yyyy-zzz, where xxxxx is the middle Julian date 
(MJD) of the observation, yyyy is the plate number, and zzz is the fiber number. 

''Abundances are calculated adopting Te=Te(0 III). 

'^Abundance are calculated adopting Te=Te(He"'"). 

''In units 10~^. 

'^In units 10"^. 



Table 4. Parameters for the best solution for He mass fraction in Table [3] 



Object'' 



Te(0 III) 



Te(He+) 



t(A3889) AHa/Ha ICF 



CGCG 007-025 (No. 1) 1.64+ 0.02 1.64+ 0.04 295 

CGCG 007-025 (No. 2) 1.65+ 0.03 1.57+ 0.04 13 

HS 0029+1748 1.28+ 0.02 1.22+ 0.03 12 

HS 0111+2115 1.11+ 0.06 1.10+ 0.03 459 

HS 0122+0743 1.76+ 0.02 1.68+ 0.05 11 



102 
51 
+96 
2 

47 
1 

+ 1 

268 
66 




0.75 

0.94 
1.16 
1.15 
1.05 



+0.16 
0.23 
+0.26 
0.59 
+0.24 
0.36 
1.54 
0.59 
+0.12 
0.50 



0.0144 0.9957 0.72E+00 

0.0001 0.9963 0.21E+01 

0.0000 0.9959 0.78E+01 

0.0427 1.0000 0.43E+00 

0.0003 0.9957 0.29E+01 



(Abridged. Table with all H II regions will be published in the online edition of ApJ) 



''The first 93 entries are HeBCD H II regions with names in alphabetical order. The remaining 271 entries are 
SDSS H 11 regions. Their names are in the format xxxxx-yyyy-zzz, where xxxxx is the middle Julian date (MJD) 
of the observation, yyyy is the plate number, and zzz is the fiber number. 



Table 5. Budget of different systematics effects in the Yp determination 



Property 




AYp 


He I emissivity 


< 


+1.7% 


re(He+) = (0.95 - 1.0)xTe(O ill) 


< 


-1.0% 


Underlying He I stellar absorption 


< 


+3.0% 


Collisional excitation of hydrogen emission lines 


< 


+1.0% 


JCF(He+ + He2+) 


< 


-1.0% 



-56- 



Table 6. Maximum Likelihood Linear Regressions 



Method 



N 



Oxygen 



Nitrogen 



Old He I cmissivitics 
Old He I emissivities 
New He i emissivities 
New He i emissivities 

Old He I emissivities 
Old He I emissivities 
New He i emissivities 
New He i emissivities 



a) 0/H=0/H[Te(0 m)], N/H=N/H[Te(0 iii)] 

93 0.2466±0.0012 + 53±9(0/H) 0.2485±0.0009 + 
364 0.2460±0.0011 + 66±7(0/H) 0.2473±0.0007 + 
93 0.2509±0.0012 + 52±9(0/H) 0.2529±0.0009 + 
364 0.2495±0.0010 + 65±7(0/H) 0.2518±0.0007 + 

b) 0/H=0/H[re(He+)], N/H=N/H[Te(He+)] 

93 0.2472±0.0012 + 43±8(0/H) 0.2489±0.0009 + 
364 0.2457±0.0010 + 56±6(0/H) 0.2475±0.0007 + 
93 0.2516±0.0011 + 40±7(O/H) 0.2532±0.0009 + 
364 0.2505±0.0010 + 50±6(O/H) 0.2522±0.0007 + 



966±157(N/H) 
1121±104(N/H) 

914±158(N/H) 
1095±106(N/H) 

824±140(N/H) 
1030± 96(N/H) 

756±135(N/H) 
958± 94(N/H) 



